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I review progress made in solving gauge theories such as collinear quantum chro- 
modynamics using light-cone Hamiltonian methods. I also show how the light- 
cone Fock expansion for hadron wavefunctions can be used to compute operator 
matrix elements such as decay amplitudes, form factors, distribution amplitudes, 
and structure functions, and how it provides a tool for exploring novel features 
of QCD. I also review commensurate scale relations, leading-twist identities which 
relate physical observables to each other, thus eliminating renormalization scale 
and scheme ambiguities in perturbative QCD predictions. 



1 Introduction 

The key challenge of nonperturbative quantum chrome-dynamics is to compute 
the spectrum of hadrons and gluonic states from first principles as well as 
determine the wavefunctions for each QCD bound state in terms of its quark 
and gluon degrees of freedom. If we had such a complete solution, then we could 
compute the quark and gluon structure functions and distribution amplitudes 
which control hard-scattering inclusive and exclusive reactions, as well as all 
of the operator matrix elements of currents which underlie electro-weak form 
factors and the weak decay amplitudes of the light and heavy hadrons. The 
knowledge of hadron wavefunctions would also provide a deep understanding 
of the physics of QCD at the amplitude level, illuminating exotic effects of the 
theory such as color transparency, intrinsic heavy quark effects, hidden color, 
diffractive processes, and the QCD van der Waals interactions. 

Solving a quantum field theory such as QCD is clearly not easy. However, 
highly non-trivial, one-space one-time relativistic quantum field theories which 
mimic many of the features of QGD have already been completely solved using 
light-cone Hamiltonian methodsu In fact, virtually any (1+1) quantum field 
theory canJpe solved using the method of Discretized Light-Cone-Quantization 
(DLCQ).qIj In DLCQ, a quantum field theory is rendered discrete in momen- 
tum space by imposing periodic or anti-periodic boundary conditions. The 
Hamiltonian H^c, which can be constructed from the Lagrangian using light- 
cone time quantization, can then be diagonalized, in analogy to Heisenberg's 
solution of the eigenvalue problem in quantum mechanics. In the one-space 
one-time theories, the diagonalization is a straightforward computational prob- 
lem, and the resulting eigenvalues and eigensolutions then provide the complete 
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spectrum of hadrons, together with their respective light-cone wavefunctions. 

A beautiful illustration of the application of light-cone quantization to 
the solution of a quantum field theory is the DLCQ analysis of "collinear" 
QCD: a variant of QCD(3 + 1) defined by dropping all of interaction terms in 
H®£ D involving transverse momentaP Even though this theory is effectively 
two-dimensional, the transversely-polarized degrees of freedom of the gluon 
field are retained as two scalar fields. Antonuccio and DalleyS have recently 
used DLCQ to solve this theory. The diagonalization of iJ|°ll mear provides 
not only the complete bound and continuum spectrum of the collinear theory, 
but it also yields the complete ensemble of light-cone Fock state wavefunctions 
needed to construct the quark and gluon structure functions of each hadronic 
and gluonic state. For example, Antonuccio and Dallcy obtain the spectrum of 
gluonia, and the polarized gluon and quark structure functions of the mesonic 
states. Although the collinear theory is a drastic approximation to physical 
QCD(3 + 1), the phenomenology of its DLCQ solutions demonstrate general 
features of gauge theory, such as the peaking of the wavefunctions at minimal 
invariant mass, color coherence, and the helicity retention of leading partons 
in the polarized structure functions at x — > 1. The DLCQ solutions of the one- 
space one-time gauge theories can be obtained for arbitrary coupling strength, 
flavors, and colors. 

The solutions to collinear QCD provide a "standard candle" or theoretical 
laboratory in which other nonperturbative methods proposed to solve QCD, 
such as lattice gauge theory, Bcthe-Salpeter methods, and various approxima- 
tions can be tested and compared. 

The fact that one actually solve a non-trivial relativistic quantum field the- 
ory in one space and one time gives hope that the full solutions to QCD(3+1) 
will eventually be accomplished. In these lectures I shall also discuss the pos- 
sibility that one can use the collinear theory as a first approximation to a 
procedure which systematically constructs the full wavefunction solutions of 
QCD(3+1). I will also outline the progress made in understanding hadrons at 
the amplitude level, using the light-cone Fock expansion as a physics tool for 
exploring novel features of QCD in hadron physics. I also review commensu- 
rate scale relations, a method which relates physical observables to each other, 
thus eliminating ambiguities due to scale and scheme ambiguities. 

2 The Light-Cone Fock Expansion 

The concept of the "number of constituents" of a relativistic bound state such 
as a hadron in quantum chromodynamics, is not only frame-dependent, but its 
value can fluctuate to an arbitrary number of quanta. Thus when a laser beam 
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crosses a proton at fixed "light-cone" time r = t + z/c = x° +x z , an interacting 
photon can encounter a state with any given number of quarks, anti-quarks, 
and gluons in flight (as long as n q — n-g = 3). The probability amplitude 
for each such n-particle state of on-mass shell quarks and gluons in a hadron 
is given by a light-cone Fock state wavefunction il) n i H (xi,k^ i ,\ i ), where the 
constituents have positive longitudinal light-cone momentum fractions 

Xi = jjl = k [ ± k * , y^ = i, (i) 

i=l 

relative transverse momentum 

n 

k ±l , ^ = 0^ ' ( 2 ) 

i=l 

and helicities A,. The ensemble {ip n /ii} of such light-cone Fock wavefunctions is 
a key concept for hadron physics, providing a conceptual basis for representing 
physical hadrons (and also nuclei) in terms of their fundamental quark and 
gluon degrees of freedomQ In the light-cone formalism, the vacuum is essentially 
trivial. Since each particle moves forward in light-cone time r with positive 
light-cone momenta fractions x% , light-cone perturbation theory is particularly 
simple and intuitive, involving many fewer diagrams than equal-time theory. 

The light-cone Fock expansion is defined in the following way: one first 
constructs the light-cone time evolution operator P~ = P Q — P z and the 
invariant mass operator Hlc = P P + — P± in light-cone gauge A + = from 
the QCD Lagrangian. The total longitudinal momentum P + = P° + P z and 
transverse momenta P± are conserved, i.e., are independent of the interactions. 
The matrix elements of Hlc on the complete orthonormal basis {|n >} of the 
free theory H LC = Hlc (<7 = 0) can then be constructed. The matrix elements 
(n | Hlc I m ) connect Fock states differing by 0, 1, or 2 quark or gluon quanta, 
and they include the instantaneous quark and gluon contributions imposed by 
eliminating dependent degrees of freedom in light-cone gauge. 

In practice it is essential to introduce an ultraviolet regulator in order to 
limit the total range of (n | Hlc \ Tn), such as a "global" cutoff in the invariant 
mass of the free Fock states: 

Ml = E < Ag lobal . (3) 

i=l 

One can also introduce a "local" cutoff to limit the change in invariant mass 
— Mm\ < AfocaJ wn ich provides spectator-independent regularization of 
the sub-divergences associated with mass and coupling renormalization. 
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The natural renormalization scheme for the coupling is ay (Q), the effective 
charge defined from the scattering of two infinitely heavy quark test charges. 
The renormalization scale can then be determined from the virtuality oLthe 
exchanged momentum, as in the BLM and commensurate scale methodsLm I 
will discuss this further in Section 18. 

In the DLCQ method , the matrix elements (n | H^q \ rnj, are made dis- 
crete in momentum space by imposing periodic or anti-periodic boundary con- 
ditions in x~ — x° — x z and x± (see Section 3). Upon diagonalization of Hlc, 
the eigenvalues provide the invariant mass of the bound states and eigenstates 
of the continuum. The projection of the hadronic eigensolutions on the free 
Fock basis define the light-cone wavefunctions. For example, for the proton, 

\P) = ^2(n\p) \n) 

n 

= ^39/p( X .£j-«>**) \ Uud ) ( 4 ) 
+^3g]/p( X i> X i) I Uud 9) H 

The light-cone formalism has the remarkable feature that the ip^ H (xi, fe_u, A c ) 
are invariant under longitudinal boosts; i.e., they are independent of the to- 
tal momentum P + , P± of the hadron. As we shall discuss below, given the 
^wj?' we can cons t ruc t an y electromagnetic or electroweak form factor from 
the diagonal overlap of the LC wavefunctions i This is illustrated in detail in 
Section 6. Similarly, the matrix elements of the currents that define quark and 
gluon structure functions can be computed from the integrated squares of the 
LC wavefunctionsEiJ 

These properties of the LC formalism are all highly-nontrivial features. In 
contrast, in equal-time formalism, the evaluation of any electromagnetic form 
factor requires the computation of non-diagonal matrix elements of bound state 
Fock wavefunctions in which the parton number can change by two units; even 
worse, matrix elements involving spontaneous pair production or annihilation 
is also required, so a complete solution of the vacuum is also required. In 
light-cone quantization, the full vacuum is also the vacuum of the free theory 
and thus is trivial. Further, the equal-time computation is only valid in one 
Lorentz frame; boosting the result to a different frame is a dynamical problem 
as complicated as solving the complete Hamiltonian problem itself. 

In general, any hadronic amplitude such as quarkonium decay, heavy 
hadron decay, or any hard exclusive hadron process can be constructed as 
the convolution of the light-cone Fock state wavefunctions with quark-gluon 
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matrix elements 

/n n / n \ / n \ 

« ti i=l i=l \ i=l / \i=l / 

x^Csi.fcxi.Ai)^ . (5) 

Here Mq A g is the underlying quark-gluon subprocess scattering amplitude, 
where the (incident or final) hadrons are replaced by quarks and gluons with 
momenta Xip + , Xip±+k±i and invariant mass above the separation scale .M 2 > 
A 2 . The LC ultraviolet regulators thus provide a LC factorization scheme for 
elastic and inelastic scattering, separating the hard dynamical contributions 
with invariant mass squared A4 2 > A 2 lobal from the soft physics with Ai 2 < 
A.giobai which is incorporated in the nonperturbative LC wavefunctions. The 
DGLAP evolution of parton distributions can be deriiied by computing the 
variation of the Fock expansion with respect to A 2 lobal l_l 

The use of the global cutoff to separate hard and soft physics is more than 
a convention; it is essential in order to correctly analyze the behavior oLdeep 
inelastic scattering structure functions in the Xbj — > 1 endpoint regimes At 
large x, the spectator constituents of the hadron target are forced to stop, 
placing the struck quark far off shell. From the stand-point of the LC Hamil- 
tonian theory, the LC energy M 2 — .M 2 becomes infinitely negative. Similarly, 
in the covariant formalism, the Feynman virtuality becomes infinitely space- 
like: k% - to 2 = - p s ) 2 - m 2 = -x(M 2 - Ml) -> -(to 2 + k%_)/(l - x) 
for x — > 1. Here x is the light-cone momentum fraction of the light quark, p s 
is the four-momentum of the remnant spectator system with p 2 = to 2 > 0. 
Thus in the large x regime, where kp becomes far off-shell, one cannot sep- 
arate the hard physics from the physics of the target wavefunction. The 
LC factorization scheme correctly isolates this phenomena. An important 
physical consequence is that DGLAP evolution is truncated; effectively, the 
starting evolution scale is of order — kp. Thus in the xbj ~ 1 regime where 
W 2 — (p + q) 2 = (1 — XBj)Q 2 /x^ is fixed, the DGLAP evolution due to the 
radiation of hard gluons is truncated and the power law behavior (1 — XBj) n 
dictated by the underlying hadron wavefunctions is maintained. It is this fact 
which allows exclusive- inclusive duality in deep inelastic lepton scattering to be 
maintained: the nominal power law (dimensional counting) behavior of struc- 
ture functions at x — > 1 and the power-law fall off of form factors at large Q 2 
match at fixed W 2 . 

The simplest, but most fundamental, characteristic of a hadron in the 
light-cone representation, is the hadronic distribution amplitude(L2l defined as 
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the integral over transverse momenta of the valence (lowest particle number) 
Fock wavcfunction; e.g., for the pion 

Mxi, Q) = J d 2 k ± ^j w (xi,k Xi ,X) (6) 

where the global cutoff A g ; f, a i is identified with the resolution Q. The distribu- 
tion amplitude controls leading-twist exclusive amplitudes at high momentum 
transfer, and it can be related to the gauge-invariant Bethe-Salpeter wavefunc- 
tion at equal light-cone time r = x + . The distribution amplitude is boost and 
gauge invariant. Its logQ evolution can be derived from the perturbatively- 
computable tail of rthp, valence light-cone wavcfunction in the high transverse 
momentum regimeta'O 

Exclusive processes are particularly challenging to compute in quantum 
chromodynamics because of their sensitivity to the unknown nonperturbative 
bound state dynamics of the hadrons. However, in some important cases, 
the leading power-law behavior of an exclusive amplitude at large momentum 
transfer can be computed rigorously in the form of a factorization theorem 
which separates the soft and hard dynamics. For example, the leading l/Q 2 
fall-off of the meson form factors can be computed as a perturbative expansion 
in the QCD coupling^M 

F M (Q 2 )= [ dx [ dy<t> M {x,Q)T H {x,y,Q 2 )<t> M {y,Q), (7) 
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where 4>m(x, Q) is the process- independent meson distribution amplitude which 
encodes the nonperturbative dynamics of the bound valence Fock state up to 
the resolution scale Q, and 

T H (x, y, Q 2 ) = ™*C F (1 + 0(a s )) (8) 

is the leading-twist perturbatively calculable subprocess amplitude 
j*q(x)q(l — x) — ► q(y)q(l — y), obtained by replacing the incident and final 
mesons by valence quarks collinear up to the resolution scale Q. The contri- 
butions from non-valence Fock states and the correction from neglecting the 
transverse momentum in the subprocess amplitude from the nonperturbative 
region are higher twist, i.e., they are power-law suppressed. The transverse 
momenta in the perturbative domain lead to the evolution of the distribution 
amplitude and to next-to- leading-order (NLO) corrections in a s . The contribu- 
tion from the endpoint regions of integration, x *~ 1 and y ~ 1, are power-law 
and Sudakov-auppressed and thus can only contribute corrections at higher 
order in 1/qN 
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The physical pion form factor must be independent of the separation scale 
Q. Again, it should be emphasized that the natural variable to make this 
separation is the light-cone energy, or equivalently the invariant mass M 2 = 

k± /x(l — x), of the off-shell partonic system@0 Any residual dependence 
on the choice of Q for the distribution amplitude will be compensated by a 
corresponding dependence of the NLO correction in Tjj- However, the NLO 
prediction for the pion form factor depends strongly on the form of the pion 
distribution amplitude as well as the choice of renormalization scale \x and 
scheme. I will discuss recent progress in eliminating such scheme and scale 
ambiguities and testing QCD in exclusive processes in Sections 18 and 19. 

A catalog of applications of light-cone Fock wavefunctions to QCD pro- 
cesses is illustrated in Fig. [j]. The light-cone expansion for a proton in terms of 
the complete set of color-singlet baryon number B = 1 free Fock states is illus- 
trated in Fig. 0a. The distribution amplitude which controls high momentum 
transfer mesonic processes is illustrated in Fig. |l|b. The meson distribution am- 
plitude 4>m(x, Q) is defined such that the invariant mass M. of the free partons 
in any intermediate state are cutoff at the ultraviolet scale Q. The relation of 
structure functions in deep inelastic lepton scattering to the integrated square 
of light-cone Fock wavefunctions is illustrated in Fig. 0c. As shown in Fig. |l|d, 
the light-cone wavefunctions provide an exact basis for the computation of the 
matrix elements of spacelike local currents in terms of diagonal overlap inte- 
grals with Xi — x[ unchanged and parton number n = n' . Figure |l|d also illus- 
trates the factorization of the nucleon form factors at high momentum transfer 
in terms of the convolution of a hard-scattering amplitude Tjj (xi, Q) (for 
the scattering the valence quarks from the initial to final direction) with the 
nucleon distribution amplitudes. Figure |I]e illustrates the application of per- 
turbative QCD factorization to the Compton scattering amplitude at large —t 
and s. Figure 0f illustrates the computation of a non- forward matrix element 
of currents; specifically virtual Compton scattering where the incident photon 
has large incident virtuality q 2 — —Q 2 . Non-diagonal n! = n — 2 Fock state 
convolutions are required when evaluating such processes in a collinear refer- 
ence frame. The computation of a weak decay matrix element B° — ► D + tu in 
terms of light-cone Fock wavefunctions is illustrated in Fig. [|g. Both diagonal 
and non-diagonal overlap integrals contribute. Finally, the perturbative QCD 
factorization of the dominant contribution to vector meson leptoproduction at 
large photon virtuality Q 2 and high energy s ^> —t, Q 2 is illustrated in Fig. 
0h. The dominant contribution arises from longitudinally polarized photons, 
and the qq pair typically has small transverse size bj_ ~ 1/Q. The coupling of 
the quark pair to the outgoing vector meson is controlled by the vector meson 
distribution amplitude. 
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Figure 1: Computation of QCD processes in terms of light-cone Fock wavefunctions. The 
various processes are described in the text. 
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3 DLCQ: A Program for Solving QCD (3+1). 

The DLCQ method i consists of diagonalizing the light-cone Hamiltonian at 
fixed x + on a free Fock basis { | n)}; i.e. the complete set of eigenstates of the 
free Hamiltonian U\ c satisfying periodic or anti-periodic boundary conditions 
in x~ . The eigenvalue problem is 

H LC | *) = M 2 | *) (9) 
(n | H LC | m) (to I *) = M 2 (n | #) (10) 

with 

kt = —m>0 P+ = —K Y>n l= K. (11) 
L L 

Here K, the "harmonic resolution," is an arbitrary positive integer. The con- 
tinuum limit corresponds to K oo. The value of length L is an irrelevant 
boost parameter in that it never appears in physics results. Since there are 
only a finite number of partitions of a given K among the positive integers 
rij with Sni = K, the number of distribution Fock states are automatically 
rendered discrete. The transverse momenta are also made discrete by choosing 
periodic or anti-periodic conditions in x±. Then 

k±i = — nj_i . (12) 

The limit on the number of states is then controlled by the global cutoff. 

The diagonalization of the light-cone Hamiltonian thus becomes the prob- 
lem of diagonalizing large Hermitian matrices, a numerical analysis problem, 
solvable by Lanczos or other methods. In the case of 1+1 dimensions, the 
problem is completely tractable, so virtually any 1+1 quantum field theory 
can be solved in this manner. 

In the case of QCD (1+1) in A + = gauge, where there are no dynamical 
gluons, the only interaction terms arise from "instantaneous" gluon exchange: 

(13) 



H\ c = 9 - 

7T 



_(k+-£+) 2 (fc++m+) 2 _ 

corresponding to i-channcl or s-channel contributions in the amplitude 

q(k + ) q(m + ) + q(£ + ) q(n + ) . (14) 

There is also a mass renormalization contribution to q(n + ) generated from 
normal ordering 

2 ™ + i 

5H= g -Y.— 2 - (15) 

m—l 
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The solution to the diagonalization of the light-cone Hamiltonian produces 
not only the mass eigenvalues of the theory, but also the eigensolutions, as 
wavefunction coefficients in the LC Fock basis: 

\^) = J2^i,k±i,Xi) |n) (16) 

n 

where each Fock component | n) has the same global and conserved quantum 
numbers as the eigenstate. The values of the light-cone momentum fractions 
are evaluated at 

,+ _*£_:*- K ~ 1 } (u) 

* p+ K [K'K'K K J [ ' 

Thus one samples the wavefunctions at rational points which approach the 
continuum theory at if — > oo. The absence of the end-points at Xi = 0, 1 
corresponds to the neglect of zero modes. Except for massless, collinear k\, 
such parton configurations are associated with infinitely massive free energy: 

Ml = ±{^) r *> (18) 

and thus exceed the global cutoff limit. Physically, the x — > limit is associated 
with partons infinitely far in rapidity from the center of mass of the bound state 
itself 

y i -Y = £nx i ; (19) 

such partons are only relevant at very large energies W 2 = (p + q) 2 in the 
computation of structure functions. On the other hand, the LC Fock wave- 
functions do not necessarily vanish at Xj = since they may correspond to soft 
gluons with m 2 and k\ = 0. In general, even the fermion distribution need 
not vanish at x — > in gauge theory since only the combination from the sum 
over states, 

"(le± -gA ± ) 2 +m 2 



E 



(20) 



has to be finite in the interacting theory. As Antonuccio and Dalley and 10 
have recently shown, the cancellation of infinities at x\ — > for fermions in 
gauge theory imposes strict "ladder relations" between Fock states with one 
or two more or less gluons in the bound state. We have also shown how this 
type of analysis leads to Regge power-law behavior of the quark distributions 
at x — > 0. 
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These interesting features are illustrated in Fig. @a and 0b which shows the 
gluon G(x) and quark Q(x) distributions of the lowest and first excited 
1 mesonic bound states of collinear QCD (N c — > oo) from DLCQB Neither 
the quark nor the gluon distributions vanish at a; — > 0. The polarization 
asymmetry for quarks and gluons with helicity aligned and anti-aligned to the 
J z spin of 1 ground state shown in Fig. |^c demonstrates the tendency 
toward helicity alignment at x — > 1. 

4 Jet Hadronization in Light-Cone QCD 

One of the goals of nonperturbative analysis in QCD is to compute jet hadroniza- 
tion from first principles. The DLCQ solutions provide a possible method to 
accomplish this. By inverting the DLCQ solutions, we can write the "bare" 
quark state of the free theory as 

ko}=£ \n)(n\q ) (21) 

where now {|n)} are the exact DLCQ eigenstates of Hlc, and (n\qo) are 
the DLCQ projections of the eigensolutions. The expansion in automatically 
infrared and ultraviolet regulated if we impose global cutoffs on the DLCQ 
basis: 

A 2 < AMI < A 2 • (22) 

where A.M 2 = .M 2 — (E.M ;) 2 . It would be interesting to study this type of jet 
hadronization at the amplitude level for the existing DLCQ solutions to QCD 
(1+1) and collinear QCD. 

5 Light-Cone Quantization and Renormalization Theory 

The renormalization procedure for LC Hamiltonian theory is well understood in 
perturbation theory. For example, mass and coupling renormalization counter 
terms can be introduced in the standard way in QED to absorb the ultraviolet 
divergences at each order of perturbation theory. An explicit method, "alter- 
nating denominators" , which provides an automatic method to construct the 
local counter terms, was usedlj to compute the lepton anomalous moments 
through order a 2 and partly through order a 3 . Lepage and I have employed an 
ultraviolet Hamiltonian renormalization scheme to derive the hard-scattering 
expansion for exclusive processes in QCD, including the evolution equations 
for the renormalized distribution amplitudes. Burkardt and LangnaiEj have 
shown that kinetic and vertex mass renormalization counter terms are needed 
to restore the full invariance structure of the theory when one uses a light-cone 
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regulation such as the global cutoff. BurkardfiS also has shown that tadpole 
diagram renormalization of gift theory is consistently handled in the LC theory 
as a zero mode component to the renormalization of the scalar particles. The 
renormalization of Yukawa theories has also been analyzed EZl Hiller, McCartor, 
and I are working on a program in which the full Lorentz-invariant structure of 
light-cone Hamiltonian theory is restored using generalized Pauli-Villars regu- 
larization. 

One of the advantages of DLCQ is that it provides a convenient infrared 
regularization of zero modes since they become discrete entities. In some model 
field theories, the zero modes take the place of the vacuum in reproducing the 
physics of spontaneous symmetry breaking. In other cases, such as the massive 
Schwinger model QED (1+1), the zero modes allow a simulation of external 
electric fields. As yet it is not clear whether LC zero roles play an essential 
role in analyzing QCD(3+1). 



6 Form Factors and Light-Cone Wavefunctions 



A critical advantage of the light-cone formalism is that the knowledge of the 
LC Fock wavefunction is sufficient to compute the elastic electroweak form 
factors. It is remarkable that all such matrix elements can be computed from 
diagonal (parton-conserving) overlap integrals of the LC Fock wavefunctions. 
In this section I will review the light-cone formalism for the .computation of 
form factors for both elementary and composite systemsli3'EJM We can choose 
light-cone coordinates with the incident lepton directed along the z direction^] 
{p ± = p°±p 3 ): 



(p + ,p ,p\) 



M 2 



-,o_ 



o, — t ; q± 



(23) 



where q 2 = —2q ■ p = —q\ and M = mi is the mass-of the composite system. 
The Dirac and Pauli form factors can be identifiecfi'E3 from the spin-conserving 
and spin- flip current matrix elements (J + = J° + J 3 ): 



= {p + q. 



M 



u 



J+(0) 



p + q, T 



p + 

J+(0) 



P ,T = 2 W) 



P, T > = -2(9! - iq 2 



F 2 (q 2 ) 



where f corresponds to positive spin projection S z 



2M ' 
i along the z axis. 



(24) 
(25) 
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Each Fock-state wave function \n) of the incident lepton is represented by 
the functions ?Pps z ( Xi ' ^i)> wnere 

k» = {k+ ,k- X) = (xp+, fc±+ , m ,JL 

\ xp + 

specifies the light-cone momentum coordinates of each constituent i = 1, . . . ,n, 
and Si specifies its spin projection S*. Momentum observation on the light 
cone requires 

n n 

^2 k±i = , ^2 x i = 1 j 

i=l i=l 

and thus < Xj < 1. The amplitude to find n (on-mass-shell) constituents 
in the lepton is then ip^> multiplied by the spinor factors us i (ki)/(k^)^ 1 ^ 2 or 
vSi (hi) / (kf) 1 / 2 for each constituent fermion or anti-fermion. The Fock state 
is off the "energy shell" : 

The quantity (fctu + m i)/ x i is t ne relativistic analog of the kinetic energy 
pf/2rrii in the Schrodinger formalism. 

The wave function for the lepton directed along the final direction p + q in 
the current matrix element is then 

wher(0 

feij = 4j + (1 - 
for the struck constituent and 

for each spectator [i ^ j). The k'± are transverse to the p + q direction with 

n 
i=l 

The interaction of the current J + (0) conserves the spin projection of the 
struck constituent fermion (u s , -f + u s )/k+ = 25 ss >. Thus from Eqs. and 
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Fi{q 2 ) = =Y / e j J[dx] [d 2 b] V^T (x,k' ± ,s) ^ (x,k ± ,s) , 



(26) 



and 



gi - iqi 
2M 



E e i \ 



d 2 fc. 



where ej is the fractional charge of each constituent. [A summation of all 
possible Fock states (n) and spins (S) is assumed.] The phase-space integration 
is 

n 

(28) 



[dx] = 8 (1 - y^ y Xj ) Y[dx, , 



and 



(29) 



Equation (^) evaluated at q 2 = with Fi(0) = 1 is equivalent to wavefunction 
normalization. The anomalous moment a = i 7 2(0)/i 7 i(0) can be determined 
from the coefficient linear in q\ — iq2 from the coefficient linear in q\ — iq2 from 
ip* +q in Eq. (§?]). In fact, 



d_ 



E 



d 



■4> 



p+q 



(30) 



(summed over spectators), we can, after integration by parts, write explicitly 

^3 



dk 



2, 



The wave function normalization is 

/ [da;] / [d 2 k±] ip^ Vyr = / [dx] / d 2 k^ij/* l ip pl = 1 



(31) 



(32) 



A sum over all contributing Fock states is assumed in Eqs. (|3l] ) and (32). We 
thus can express the anomalous moment in terms of a local matrix element at 



1G 



zero momentum transfer. It should be emphasized that Eq. (|3l]) is exact; it is 
valid for the anomalous element of any spin-i system. nr _. 

In the case of the electron's anomalous moment to order a in QEDptJ 
the contributing intermediate Fock states are the electron-photon states with 
spins — 5,1) and 5,— l): 

elJx { V2^^ (1-1)^1-1,1)) 

^ e Jxi — y 27 12/7 (33) 

and 

a,. J-A'^lH-Wi)) (34 ) 



The quantities to the left of the curly bracket in Eqs. ( |33| ) and (|34|) are the 
matrix elements of 

— u u 

■7-e 7-ZTT75 and /_4.M/ g 7-e 



(p+ - fc+) 1 /2 ' (p+)Va (p+)!/2 ' (p+ _ fc+)l/2 ' 

respectively, where e = e^p = ±(l/v / 2)(5; ± zy), e_;ii 1 = 0, e+ = in the 
light-cone gauge for vector spin projection S z = ±lli§Ea For the sake of gen- 
erality, we let the intermediate lepton and vector boson have mass fh and A, 
respectively. 

Substituting @ and Q into Eq. @, one finds that only the |-|,l) 
intermediate state actually contributes to a, since terms .which involve differ- 
entiation of the denominator of ippi cancel. We thus havla 

a = 4Me* [ **± f to (35 ) 
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[M 2 - {k]_ + m 2 )/(l -x)- (k 2 ± + A 2 )/a 



or 



r a f 1 M[m-M(l-x)]x(l-x) 

-I m 2 x + \ 2 (l-x)-M 2 x{l-x) ' 1 ' 



7T 



(I 



which, in the case of QED (m = M, A = 0) gives the Schwinger results a 
a/2-K. 

The general result (|3l]) can also be written in matrix form: 



^ = - E e ^ / J ^ • ^ , (37) 
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where S is the spin operator for, the total system and Lj_ is the generator of 
"Galiean" transverse boosts HL3 on the light cone, i.e., S±_ ■ L± = (5 + L_ + 



'HrY 

S-L + )/2 where S± = (Si ± iS?) is the spin-ladder operator and 



*±-i>(sb*'sb) < 38 > 

(summed over spectators) in the analog of the angular momentum operator 
p* x r*. Equation ( pl| ) can also be written simply as an expectation value in 
impact space. 

The results given in Eqs. (p6|), (p7|), and ( |3l| ) are also valid for calculating 
the anomalous moments and form factors of hadrons in quantum chromody- 
namics directly from the quark and gluon wave functions tp(k±, x, S). These 
wave functions can also be used to construct the structure functions and dis- 
tribution amplitudes-jszhich control large momentum transfer inclusive and 
exclusive processes. The charge radius of a composite system can also be 
written in the form of a local, forward matrix clement: 



dFi(g 2 
dq 2 




ej J [dx] [d 2 k ± ] r pA VfcT ■ (39) 



7 Magnetic and Electroweak Moments of Nucleons in the Light- 
Cone Formalism 



The use of covariant kinematics leads to a number of striking conclusions for 
the electromagnetic and weak moments of nucleons and nuclei. For exam- 
ple, magnetic moments cannot be written as the naive sum 7^ = S 01 
the magnetic moments of the constituents, except in the nonrelativistic limit 
where the radius of the bound state is much larger than its Compton scale: 
RaMa 3> 1. The deuteron quadrupole moment is in general nonzero even if 
the nucleon-nucleon bound state has no -D-wave componentEj The breakdown 
of simple additivity for moments and the contradictions with the traditional 
nonrelativistic formalism, even for weak binding, is due to the fact that the 
so-called "static" moments must be computed as transitions between states of 
different momentum and + , with q^ — > 0. Thus one must construct 
current matrix elements between boosted states. The Wigner boost-generates 
nontrivial corrections to the current interactions of bound systems. c3 Remark- 
ably, in the case of the deuteron, both the quadrupole and magnetic moments 
become equal to that of the Standard Model in the limit MdRd ^ 0. In this 
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limit, the three form factors of the deuteron have the same ratios as do those 
of the W boson in the Standard ModelcJ 

One can also use light-cone methods to show that the proton's magnetic 
moment \x p and its axial- vector coupling qa have-a relationship independent of 
the specific form of the light-cone wavefunctionc3 At the physical value of the 
proton radius computed from the slope of the Dirac form factor, Ri = 0.76 fm, 
one obtains the experimental values for both fi p and qa ; the helicity carried 
by the valence u and d quarks are each reduced by a factor ~ 0.75 relative to 
their nonrelativistic values. At infinitely small radius R p M p — > 0, fi p becomes 
equal .to. the Dirac moment, as demanded by the Drell-Hearn-Gerasimov sum 
ruleEZrEa Another surprising fact is that as R\ — > the constituent quark 
hclicities become completely disoriented and qa — * 0. 

One can understand the origins of the above universal features even in 
an effective three-quark light-cone Fock description of the nucleon. In such a 
model, one assumes that additional degrees of freedom_(including zero modes) 
can be parameterized through an effective potentiallia After truncation, one 
could in principle obtain the mass M and light-cone wavefunction of the three- 
quark bound-states by solving the Hamiltonian eigenvalue problem. It is rea- 
sonable to assume that adding .more quark and gluonic excitations will only 
refine this initial approximationEj In such a theory the constituent quarks will 
also acquire effective masses and form factors. 

Since we do not have an explicit representation for the effective potential in 
the light-cone Hamiltonian P~ s for three quarks, we shall proceed by making 
an Ansatz for the momentum-space structure of the wavefunction Even 
without explicit solutions of the Hamiltonian eigenvalue problem, one knows 
that the helicity and flavor structure of the baryon eigenfunctions will reflect 
the assumed global SU(6) symmetry and Lorentz invariance of the theory. 
As we will show below, for a given size of the proton the predictions and 
interrelations between observables at Q 2 = 0, such as the proton magnetic 
moment fi p and its axial coupling-.^ , turn out to be essentially independent 
of the shape of the wavefunctioncB 

The light-cone model given by Ma0 and by SchlumpfEl provides a frame- 
work for representing the general structure of the effective three-quark wave- 
functions for baryons. The wavefunction \P is constructed as the product of a 
momentum wavefunction, which is spherically symmetric and invariant under 
permutations, and a spin-isospin wave function, which is uniquely determined 
by SU(6)-symmetry requirements. A Wigner-Melosh rotationoo is applied to 
the spinors, so that-the, wavefunction of the proton is an eigenfunction of J and 
J z in its rest frameEJoESI To represent the range of uncertainty in the possible 
form of the momentum wavefunction, one can choose two simple functions of 
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the invariant mass M of the quarks: 



^u.o.{M 2 ) = Nn.o.eM~M 2 /2p 2 ), (40) 
^Pow Cr (M 2 ) = N Powci (l+M 2 /p 2 y p , (41) 

where (3 sets the characteristic internal momentum scale. Perturbative QGD 
predicts a nominal power-law fall off at large k± corresponding to p = 3.5lj 
The Melosh rotation insures that the nucleon has j = \ in its rest system. It 
has the matrix representation Ej 

m + XjM -ilf -(nx kj) 

R ^> k ^ = ^ m + Xl MY + k it (42) 

with n = (0,0, 1), and it becomes the unit matrix if the quarks are collinear, 
RM(xi,0,m) = 1. Thus the internal transverse momentum. dependence of the 
light-cone wavefunctions also affects its helicity structureEa 

As we showed in Section 6, the Dirac and Pauli form factors F\(Q 2 ) and 
F2 (Q 2 ) of the nucleons are given by the spin-conserving and the spin- flip matrix 
elements of the vector current Jy (at Q 2 = — q 2 )a 

Fi(Q 2 ) = (p + 9,T|4b,T), (43) 
(Qi-*Q2)F 2 (Q 2 ) = -2M(j, + q,U4\p,l) ■ (44) 

We then can calculate the anomalous magnetic moment a = limQ2^ ^2(Q 2 )-[| 
The same parameters as given by Schlumpf c3 are chosen, namely m — 0.263 
GeV (0.26 GeV) for the up (down) quark masses, = 0.607 GeV (0.55 GeV) 
for V'Power (V'h.o.), and p — 3.5. The quark currents are taken as elementary 
currents with Dirac moments tt 3 -. All of the baryon moments are well- fit 
if one takes the strange quark mass as 0.38 GeV. With the above values, the 
proton magnetic moment is 2.81 nuclear magnetons, and the neutron magnetic 
moment is —1.66 nuclear magnetons. (The neutron value can be improved by 
relaxing the assumption of isospin symmetry.) The radius of the proton is 0.76 
fm, i.e., MpRi — 3.63. 

In Fig. |^(a) we show the functional relationship between the anomalous 
moment a p and its Dirac radius predicted by the three-quark light-cone model. 
The value of 



?2 _ MQ 2 ) 



R{ = -6- 



(45) 

2 =o 



dQ 2 

is varied by changing (3 in the light-cone wavefunction while keeping the quark 
mass m fixed. The prediction for the power-law wavefunction V'Power is given 

*The total proton magnetic moment is fi p = ^tt (1 + &p)- 
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Figure 3: (a). The anomalous magnetic moment of the proton a p = F2(0) as a function 
of its Dirac radius M P R\ in Compton units, (b). The axial vector coupling of the neutron 
to proton beta-decay as a function of M p Ri. In each figure, the broken line is computed 
from a wavefunction with power-law fall off and the solid curve is computed from a Gaussian 
wavefunction. The experimental values at the physical proton Dirac radius are indicated by 

the dotted line. 
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by the broken line; the continuous line represents i/'h.o.- Figure ^(a) shows 
that when one plots the dimensionless observable a p against the dimension- 
less observable MRi the prediction is essentially independent of the assumed 
power-law or Gaussian form of the three-quark light-cone wavefunction. Dif- 
ferent values of p > 2 also do not affect the functional dependence of a p (M p Ri) 
shown in Fig. ||(a). In this sense the predictions of the three-quark light-cone 
model relating the Q 2 — > observables are essentially model-independent. The 
only parameter controlling the relation between the dimensionless observables 
in the light-cone three-quark model is m/M p which is set to 0.28. For the phys- 
ical proton radius M p R\ = 3.63 one obtains the empirical value for a p = 1.79 
(indicated by the dotted lines in Fig. ||(a)). 

The prediction for the anomalous moment a can be written analytically 
as a = (7y)a NR , where a NR = 2M p /3m is the nonrelativistic (R — > oo) value 
and 7y is given asE3 



3m 
M 



(1 - x 3 )M(m + x 3 M) - fcj 3 /2 
(m + x 3 M) 2 + k 2 ±3 



(46) 



The expectation value (7^) is evaluated asfl 

I[d 3 khv\^\ 2 



(7v 



(47) 



Let us now take a closer look at the two limits R — > 00 and R — > 0. 
In the nonrelativistic limit we let (3 — > and keep the quark mass m and 
the proton mass M p fixed. In this limit the proton radius Ri — > 00 and 
a p — > 2M p /3m — 2.38, since (7^) — * l|] Thus the physical value of the 
anomalous magnetic moment at the empirical proton radius M p R\ — 3.63 is 
reduced by 25% from its nonrelativistic value due to relativistic recoil and 
nonzero fc^.| 

To obtain the ultra-relativistic limit we let f3 — > 00 while keeping m fixed. 
In this limit the proton becomes pointlike, M p R\ — > 0, and the internal trans- 
verse momenta k± — > 00. The anomalous magnetic momentum of the proton 

tHere [d 3 k] = dk\dk2dkzi(k\ + ki + /C3). The third component of k is defined as k^i = 

^(xiM x-M~* )' This measure differs from the usual one usedlij by the Jacobian FJ d $ l 

which can be absorbed into the wavefunction. 

tThis differs slightly from the usual nonrelativistic formula 1 + a = £2. Me due ^ Q the 

J Z Jq e TUg 

nonvanishing binding energy which results in M p ^ 3m q . 

§The nonrelativistic value of the neutron magnetic moment is reduced by 31%. 
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goes linearly to zero as a ^JLA3M p Ri since (tv) — > 0. Indeed, the Drcll- 
Hearn-Gerasimov sum ruleEZEa demands that the proton magnetic moment 
become equal to the Dirac moment at small radius. For a spin-i system 



a M 2 f°° ds 

a 



2n 2 a 



f -[a P (s)-a A (3)] , (48) 

Js th S 



where (Jp( A \ is the total photo absorption cross section with parallel (anti- 
parallel) photon and target spins. If we take the point-like limit, such that the 
threshold for inelastic excitation becomes infinite while the mass of the system 
is kepLfinite, the integral over the photo absorption cross section vanishes and 
a = OlJ In contrast, the anomalous magnetic moment of the proton does not 
vanish in the nonrelativistic quark model as R — > 0. The nonrelativistic quark 
model does not reflect the fact that the magnetic moment of a baryon is derived 
from lepton scattering at nonzero momentum transfer, i.e., the calculation of 
a magnetic moment requires knowledge of the boosted wavefunction. The 
Melosh transformation is also essential foLderiving the DHG sum rule and 
low-energy theorems of composite systemsEd 

A similar analysis can be performed for the axial- vector coupling measured 
in neutron decay. The coupling gA is given by the spin-conserving axial current 



J_4 matrix element 



g A {Q) = (j>A \JX\pA) ■ (49) 



The value for gA can be written as gA — , with being the nonrel- 

ativistic value of gA and with ja given byEZIH 



(to + x 3 M) 2 — k 2 
(to + x 3 M) 2 + k ±3 



1 A (xi, k±i, to) = - ■ 2 . (50) 



In Fig. 0(b) the axial-vector coupling is plotted against the proton radius 
M p R\. The same parameters and the same line representation as in Fig. [|(a) 
are used. The functional dependence of g a(M p Ri) is also found to be indepen- 
dent of the assumed wavefunction. At the physical proton radius M p R\ — 3.63, 
one predicts the value gA — 1-25 (indicated by the dotted lines in_Fig. ||(b)), 
since (j A ) = 0.75. The measured value is g A = 1.2573 ± 0.0028E3 This is a 
25% reduction compared to the nonrelativistic SU(6) value 17,4 = 5/3, which 
is only valid for a proton with large radius R\ ^_1/M p . The Melosh rotation 
generated by the internal transverse momentum^ spoils the usual identifica- 
tion of the 7 + 75 quark current matrix element with the total rest-frame spin 
projection s z , thus resulting in a reduction of gA- 

Thus, given the empirical values for the proton's anomalous moment a p 
and radius M p Ri, its axial- vector coupling is automatically fixed at the value 
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gA = 1.25. This is an essentially model-independent prediction of the three- 
quark structure of the proton in QCD. The Melosh rotation of the light-cone 
wavefunction is crucial for reducing the value of the axial coupling from its 
nonrelativistic value 5/3 to its empirical value. The near equality of the ratios 
9aI 9a{R\ —> oo) and a p /a p (Ri —* oo) as a function of the proton radius 
i?i shows the wave-function independence of these quantities. We emphasize 
that at small proton radius the light-cone model predicts not only a vanishing 
anomalous moment but also lim^^o gA(M p Ri) = 0. One can understand this 
physically: in the zero radius limit the internal transverse momenta become 
infinite and the quark helicities become completely disoriented. This is in 
contradiction with chiral models, which suggest that for a zero radius composite 
baryon one should obtain the chiral symmetry result gA = 1. 

The helicity measures Am and Ad of the nucleon each experience the same 
reduction as does gA due to the Melosh effect. Indeed, the quantity Ag is 
defined by the axial current matrix element 

Ag = {p, T |<Z7 + 75<?b, T> i (51) 

and the value for Aq can be written analytically as Aq — (^A)Aq NB -, with 
Ag NR being the nonrelativistic or naive value of Ag and ja given by Eq. ([sof). 

The light-cone model also predicts that the quark helicity sum AE = 
Am + Ad vanishes as the proton radius R\ becomes small. Note that AS 
depends on the proton size, and it should not be identified as the vector sum 
of the rest-frame constituent spins. The rest-frame spin sum is not a Lorentz 
invariant for a composite systemESI Empirically, one can measure Aq from the 
first moment of the leading-twist polarized structure function gi(x,Q). In the 
light-cone and parton model descriptions, Ag = J Q dx[q^ (x) — q^(x)], where 
q^ (x) and g^(a;) can be interpreted as the probability for finding a quark or 
antiquark with longitudinal momentum fraction x and polarization parallel or 
anti-parallel to the proton helicity in the proton's infinite momentum frameEj 
[In the infinite momentum frame there is no distinction between the quark 
helicity and its spin projection s z .) Thus Ag refers to the difference of helicities 
at fixed light-cone time or at infinite momentum; it cannot be identified with 
q{s z = +i) — q(s z = — |), the spin carried by each quark flavor in the proton 
rest frame in the equal-time formalism. 

Thus the usual SU(6) values Au NR = 4/3 and Ad NR = -1/3 are only 
valid predictions for the proton at large MR\. At the physical radius the 
quark helicities are reduced by the same ratio 0.75 as is gA/g™ due to the 
Melosh rotation. Qualitative arguments for such a reduction have been given 
elsewhereeM For M P R\ = 3.63, the three-quark model predicts Au = 1, 
Ad = -1/4, and AS = Am + Ad = 0.75. Although the gluon contribution 
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AG = in our model, the general sum r ule 



^AE + AG + L z = i (52) 

is still satisfied, since the Melosh transformation effectively contributes to L z . 

Suppose one adds polarized gluons to the three-quark light-cone model. 
Then the flavor-singlet quark-loop radiative corrections to the gluon propaga- 
tor wilLgive an anomalous contribution 8(Aq) = — t^-AG to each light quark 
hclicitya The predicted value of gA = Ait — Ad is of course unchanged. For 
illustration we shall choose fjAG = 0.f5. The gluon-enhanced quark model 
then gives values which agree well with the present experimental values. 

In summary, one sees that relativistic effects are crucial for understanding 
the spin structure of nuclcons. By plotting dimcnsionless observables against 
dimensionless observables, we obtain relations that are independent of the 
momentum-space form of the three-quark light-cone wavefunctions. For exam- 
ple, the value of gA — I -25 is correctly predicted from the empirical value of the 
proton's anomalous moment. For the physical proton radius M P R\ = 3.63, the 
inclusion of the Wigner-Melosh rotation due to the finite relative transverse 
momenta of the three quarks results in a ~ 25% reduction of the nonrelativis- 
tic predictions for the anomalous magnetic moment, the axial vector coupling, 
and the quark helicity content of the proton. At zero radius, the quark he- 
licities become completely disoriented because of the large internal momenta, 
resulting in the vanishing of gA and the total quark helicity AS. 

8 Constructing Hadron Wavefunctions in Light-Cone Quantized 
QCD 

Our ultimate goal is to actually calculate the light cone wavefunctions of the 
hadrons. In the next two sections, I will discuss possible methods in which one 
can obtain constraints and determine important properties of the wavefunc- 
tions, even in the absence of explicit solutions. 

A remarkable feature of collinear QCD is that although the theory is ef- 
fectively one-space and one-time, one still retains the two physical degrees of 
freedom from the transversely-polarized gluons. Thus the spectrum of collinear 
QCD contains gluonium states, as well as gluonic quanta in the higher Fock 
states of the mesons and baryons eigenstates of the theory. We have also seen 
that some of the features of the structure functions of hadrons in collinear 
QCD match well to the phenomenological features of QCD [3+ 1] such as the 
helicity retention of the leading constituents at large x — > 1 in the polarized 
structure functions. 
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Recently Antonuccio, Pinsky and I have investigated the possibility that 
one may be able to construct useful models of the light-cone wavcfunctions 
of QCD[3+1] by extension of the collinear QCD solutions. We have been 
considering two methods: 

(1) Minimal Subtraction. Let us ignore the complications of spin and write 
the solution to the n-quark/antiquark light-cone wavefunction of a hadron in 
the collinear theory in the form 



1 r 



9 2 Al 



where the functional dependence of the operator ^f n in the field variable A ± 
connects Fock states of different gluon number. The natural generalization of 
this dependence to the transverse space dependence is 



mf + (k _l 



It is interesting to note that the mechanical light-cone kinetic energy 

mf + {V \ - gl \) 2 



is the essential variable which controls the dynamics of gauge theory. This 
is in agreement with the fact that in laser physics the effective mass of an 
electron in an intense laser beam is "^effective = m l + e 2 A 2 . The laser analog 
also suggest that a classical approximation to the gauge field may be useful 
when the particle number is high. This could be appropriate when analyzing 
the physics of small x. 

(2) The Light- Cone Lippmann-Schwinger Equation. In principle, we can 
also construct the wavefunctions of QCD(3+1) starting with collinear QCD(1+1) 
solutions by systematic perturbation theory in AH, where AH contains the 
terms linear and quadratic in the transverse momenta k ±i which are neglected 
in the Hamilton H of collinear QCD. We can write the exact eigensolution of 
the full Hamiltonian as 



3+1) 



where 



1 



V'u+i) 



l 



l 



M 2 -H + ie 



Aff^i+i) 



M 2 -H + ie M 2 -H + 



+ 



AH- 



1 



M 2 -H + ie M 2 -H { 
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can be represented as the continued iteration of the Lippmann Schwinger re- 
solvant. Note that the matrix (M 2 — Ho) -1 is known to any desired precision 
from the DLCQ solution of collinear QCD. 

In each of these methods, the resulting wavefunction can be considered 
as approximate solution to 3+1 QCD hadron wavefunctions which could be 
subsequently improved by variational or other methods. 



9 Determining the Far Off-Shell Behavior of Hadron LC Wavefunc- 
tions 



In many cases of physical interest we are specifically interested in the behav- 
ior of the LC wavefunctions in the far-off-shell domain where M? n exceeds the 
global cutoff. This occurs for large parton momenta, k\ — ► oo, Xi — > 1 and for 
massive quanti-antiquark fluctuations. In fact, in such domains, we can con- 
struct the wavefunction pcrturbatiyely and obtain rigorous QCD predictions. 
The basic method is as followst^ Suppose we can solve 



M z 



(/■<} 



(53) 



in the soft-domain with M? n < y 2 , where /i 2 is of the order of a few GeV 2 



Then 



C/0\ = V" 



/ 



|n) 9(ji 2 -M^) = P^ \ $) 



(54) 



We can also define the complement projection operator with +Q^' — 
1. The full solution satisfies H | tp) = M 2 | ip) or 

PHP | i/j) + PHQ | V) = M 2 P | tp) 



Therefore 



QHP | i>) + QHQ | V) = M 2 Q \ 
1 



IV>) 



1> 



Cm) 



QH ± 



1> 



Cm) 



(55) 



(56) 



M 2 - QHQ 

where only high mass (.M 2 > /i 2 ) perturbatively calculable intermediate states 
appear. An example is the behavior of the valence wavefunction at large in- 
ternal transverse momentum. One finds at large k±_ = qj_ 

il>(x,qx) = \r— T [ dy [ d 2 e ± V(x,q ± ;y,i±)^ ) (y,e±) 



x(X— x) 



si 



a s (q±) / dyVi g (x,y) ip(y,q ± ) 



(57) 
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where 



rQ2 d 3 e, 

"Ul-Q") = I ^ 1>(2)(v,t±) (58) 



is the meson distribution amplitude, which takes the role of the wavefunction 
"at-the-origin" in analogous nonrelativistic calculations. The simple a s {q 2 L ) / 1 q\ 
fall-off of the valence wavefunction is the key input to the derivation of Himpn- 
sional counting rules for form factors and other exclusive processes in QCDEIf 
One can also derive the evolution equation for 

aiogQ 2 1 ' 

from this result and the simple properties of the one-gluon exchange kernel. 



10 Structure Functions at the End-Point 

The x — > 1 behavior of quark and gluon wavefunctions is controlled by far- 
off-shell configuration and thus can be analyzed perturbatively. The dominant 
contributions come from the lowest Fock states which contain the partons. For 
example the quark distribution in the nucleon can be computed from two iter- 
ations of the gluon exchange kernel which transfers the light-cone momentum 
to the struck quark from the two spectators which are required to stop. The 
result is a nominal power-law fall-off at x — > 1: q-[(x) *~ (1 — a;) 3 for quarks 
with helicity aligned with the proton and q± ~ (1 — x) 5 for quarks anti-aligned. 
Similarly the | qqqg) Fock state yields g-\{x) ~ (1 — x) 4 for gluons with helicity 
aligned and g± ~ (1 — a;) 6 for gluons helicity anti-aligned. Thus the partons 
with x — > 1 tend to have the same sign helicity as the bound stateKl 

In the absence of full solutions to the light-cone wavefunctions, one can 
construct a simple phenomenology of the polarized and anti-polarized structure 
functions by imposing a smooth connection between the perturbative QCD 
constraints at x — ► 1, and Regge behavior at x — ► 0, and the momentum 
and Bjorken sum rules. A complete discussion can be found in the literatures 
Recently Leader et al. have shown that this parametrization agrees remarkably 
well with the available data from SLAC and CERNej Figure [| shows the 
nominal form of the helicity distributions Aq(x) = q-\(x) — qi(x), for the valance 
quarks and gluons in the proton. 
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11 Intrinsic Hardness 

The light-cone wavefunctions contain high fluctuations of arbitrary mass; i.e. 
nonzero probabilities for massive pairs, massive sea quarks, etc. These fluc- 
tuations are of two types: extrinsic qg, gg or qq high mass pairs which are 
associated with the substructure of the constituents and are contained in or- 
dinary DGLAP evolution; and intrinsic functions which are due to the physics 
of the bound state wavefunction itself. For example intrinsic sea quark pairs 
QQ arise from diagrams which are interconnected to the valence quarks of the 
bound state hadron and thus depend on the valence quark correlations. It is 
easy to see that the leading perturbative contribution to intrinsic pairs falls as 

* mT~ (60) 

where M is the pair mass and T is the leading anomalous dimension associated 
with the valence wavefunction. The probability of any configurations with 
M 2 > Ml is then 

a 4+2r 

^ 2 >^)~^r, (61) 
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which implies a remarkably slow fall-off for large off-shell fluctuations. The 
result is universal for any intrinsic parton pair 



M 2 + P 2 




+ 



m\ + k 



'±2 



(62) 



Xl + X 2 



with 



k±i = xiP± + k± k±2 = x 2 P± - k± ; 



(63) 



i.e. for large relative k± and/or large quark mass. Hoyer and lc3 call this 
"intrinsic hardness." In the case of intrinsic charm or bottom pairs in the 
nucleoli, the LC wavefunction is maximized when A4^ is minimized; i.e. for 



urations occur at equal rapidity; i.e. where the heavy partons have highest 
momentum fractions. This is in contrast to the usual extrinsic sea quark which 
are subconstituents of the gluons and have low x. The extrinsic quarks evolve 
rapidly with a probability increasing as power of log Q 2 j M 2 . 

It is thus important to distinguish two types of quark and gluon contribu- 
tions to the nucleon sea measured in deep inelastic lepton-nucleon scattering: 
"extrinsic" and "intrinsic" . e3 The extrinsic sea quarks and gluons are created 
as part of the lepton-scattering interaction and thus exist over a very short 
time At ~ 1/Q. These factorizable contributions can be systematically de- 
rived from the QCD hard bremsstrahhmg and pair-production (gluon-splitting) 
subprocesses characteristic of leading twist perturbative QCD evolution. In 
contrast, the intrinsic sea quarks and gluons are multi-connected to the valence 
quarks and exist over a relatively long lifetime within the nucleon bound state. 
Thus the intrinsic qq pairs can arrange themselves together with the valence 
quarks of the target nucleon into the most energetically-favored meson-baryon 
fluctuations. 

Another interesting distinction between extrinsic and intrinsic sea quarks 
is that due to nonperturbative effects, the intrinsic contributions are generally 
not symmetric for sea quark and anti-quarks. For example, in the muonium 
atom (/x + e _ ) an intrinsic t + t~ pair would be asymmetric since the r + tends to 
be attracted to the electron and the t~ tends to be attracted to the opposite- 
sign muon. Thus the t~ would be expected to have a higher (x) than the r + . 
It is also possible to consider the nucleon wavefunction at low resolution as a 
fluctuating system coupling to intermediate hadronic Fock states such as non- 
interacting meson-baryon pairs. The most important fluctuations are most 



mi. 



(64) 



Xj_i = 




Thus the maximal intrinsic charm-bottom config- 
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likely to be those closest to the energy shell and thus have minimal invariant 
mass. For example, the coupling of a proton to a virtual K + A pair pro- 
vides a specific source of intrinsic strange quarks and antiquarks in the proton. 
Since the s and s quarks appear in different configurations in the lowest-lying 
hadronic pair states, their helicity and momentum distributions are distinct. 
Ma and I C3 have used an intermediate meson-baryon fluctuation model to 
model the possible s versus s and c versus c asymmetries of the intrinsic dis- 
tributions of the nucleon. We utilize a boost-invariant light-cone Fock state 
description of the hadron wavefunction which emphasizes multi-parton con- 
figurations of minimal invariant mass. We find that such fluctuations predict 
a striking sea quark/antiquark asymmetry in the corresponding momentum 
and helicity distributions in the nucleon structure functions. In particular, the 
strange and anti-strange distributions in the nucleon generally have completely 
different momentum and spin characteristics. The helicity structure of the in- 
trinsic ss is strongly asymmetric: the s quark from a A(uds) K(us) is aligned 
with the A helicity and (because of parity) is 100% anti-aligned with the nu- 
cleon spin. On the other hand, the s from the pseudoscalar kaon is unaligned. 
Ma and I have shown that this picture of quark and antiquark asymmetry in 
the momentum and helicity distributions of the nucleon sea quarks has support 
from a number of experimental observations, and we have suggested processes 
to test and measure this quark and antiquark asymmetry in the nucleon sea. 

11.1 Phenomenological Consequences of Intrinsic Charm and Bottom 

Microscopically, the intrinsic heavy-quark Fock component in the 7r~ wave- 
function, \udQQ), is generated by virtual interactions such as gg — > QQ where 
the gluons couple to two or more projectile valence quarks. The probability 
for QQ fluctuations to exist in a_Light hadron thus scales as a 2 s (mQ) / nig rela- 
tive to leading-twist productionCJ This contribution is therefore higher twist, 
and power-law suppressed compared to sea quark contributions generated by 
gluon splitting. When the projectile scatters in the target, the coherence of 
the Fock components is broken and its fluctuations can hadronize, forming 
new hadronic systems from the fluctuations l2I For example, intrinsic cc fluctu- 
ations can be liberated provided the system is probed during the characteristic 
time At = 2pi a b/M?- that such fluctuations exist. For soft interactions at mo- 
mentum scale 11, the intrinsic heavy .quark cross section is suppressed by an 
additional resolving factor cx /i 2 /mgtil The nuclear dependence arising from 
the manifestation of intrinsic charm is expected to be a a ~ oatA 2 / 3 , charac- 
teristic of soft interactions. 

In general, the dominant Fock state configurations are not far off shell 
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and thus have minimal invariant mass M 2 = '^2 i m±/xi where m± is the 
transverse mass of the i th particle in the configuration. Intrinsic QQ Fock 
components with minimum invariant mass correspond to configurations with 
equal-rapidity constituents. Thus, unlike sea quarks generated from a single 
parton, intrinsic heavy quarks tend_.to carry a larger fraction of the parent 
momentum than do the light quarksEJ In fact, if the intrinsic QQ pair coalesces 
into a quarkonium state, the momentum of the two heavy quarks is combined 
so that the quarkonium state will carry a significant fraction of the projectile 
momentum. 

There is substantial evidence for the existence of intrinsic cc fluctuations 
in the wavefunctions of light hadrons. For example, the charm structure func- 
tion of the proton measured by EMCia significantly larger than that predicted 
by photon-gluon fusion at large XBjf^ Leading charm production in ttN, 
hyperon-7V collisions also requires a charm source beyond leading twist t 
The NA3 experiment has also shown that the single Jj ip cross section at large 
xf is greater than expected from gg and qq production^ The nuclear de- 
pendence of this forward component is diffractive-like, as expected from the 
BHMT mechanism. In addition, intrinsic charm may account for the anoma- 
lous longitudinal polarization of the J/tp at large xf seen in irN — > J/i/jX 
interactionsE3 Further theoretical work is needed to establish that the data on 
direct J r /ip and xirPfoduction can be described using the higher-twist intrinsic 
charm mechanismEJ . . 

A recent analysis by Harris, Smith and Vogt of the excessively large 
charm structure function of the proton at large x as measured by the EMC 
collaboration at CERN yields an estimate that the probability P C c that the 
proton contains intrinsic charm Fock states is of the order of 0.6% ± 0.3%. In 
the case of intrinsic bottom, perturbative QCD scaling predicts 



more than an order of magnitude smaller. We can speculate that if super- 
partners of the quarks or gluons exist they must also appear in higher Fock 
states of the proton, such as \uud gluino gluino). At sufficiently high energies, 
the diffractive excitation of the proton will produce these intrinsic quarks and 
gluinos in the proton fragmentation region. Such supersymmetric particles can 
bind with the valence quarks to produce highly unusual color-singlet hybrid 
supersymmetric states such as \uud gluino) at high xp- The probability that the 
proton contains intrinsic gluinos or squarks scales with the appropriate color 
factor and scales inversely with the heavy particle mass squared relative to the 
intrinsic charm and bottom probabilities. This probability is directly reflected 
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in the production rate when the hadron is probed at a hard scale Q which is 
large compared to the virtual mass M. of the Fock state. At low virtualitics, 
the rate is suppressed by an extra resolution factor of Q 2 / M 2 . The forward 
proton fragmentation regime is a challenge to instrument at HERA, but it may 
be feasible to tag special channels involving neutral hadrons or muons. In the 
case of the gas jet fixed-target ep collisions such as at HERMES, the target 
fragments emerge at low velocity and large backward angles, and thus may be 
accessible to precise measurement. 

Double Quarkonium Hadroproduction It is quite rare for two charmo- 
nium states to be produced in the same hadronic collision. However, the NA3 
collaboration has measured a double J/ip production rate significantly above 
background in multi-muon events with 7r _ beams aL laboratory momentum 
150 and 280 GeV/c and a 400 GeV/c proton beam. B The relative double to 
single rate, / , is (3 ± 1) x 10~ 4 for pion-induced production, where er^ 
is the integrated single ip production cross section. A particularly surprising 
feature of the NA3 ir~N — ► ipipX events is that the laboratory fraction of the 
projectile momentum carried by the ipip pair is always very large, > 0.6 
at 150 GeV/c and > 0.4 at 280 GeV/c. In some events, nearly all of the 
projectile momentum is carried by the ipip system! In contrast, perturbative 
gg and qq fusion processes are expected to produce central ipip pairs, centered 
around the mean value, (x^) w 0.4-0.5, in the laboratory. The predicted 
ipip pair distributions from the intrinsic charm model provide a natural expla- 
nation of the strong forward production of double J / ip hadroproduction, and 
thus gives strong phenomcnological support for the presence of intrinsic heavy 
quark states in hadrons. 

It is clearly important for the double J/ip measurements to be repeated 
with higher statistics and at higher energies. The same intrinsic Fock states 
will also lead to the production of multi-charmed baryons in the proton frag- 
mentation region. The intrinsic heavy quark model can also be used to predict 
the features of heavier quarkonium hadroproduction, such as TT, Tip, and 
(cb) (cb) pairs. Predictions for these events have been given by Ramona Vogt 
and myself. 

Leading-Particle Effect in Open Charm Production According to PQCD 
factorization, the fragmentation of a heavy quark jet is independent of the pro- 
duction process. However, there are strong correlations between the quantum 
numbers of D mesons and the charge of the incident pion beam in ttN — > DX 
reactions. This effect can be explained as being due to the coalescence of the 
produced intrinsic charm quark with co-moving valence quarks. The same 
higher-twist recombination effect can also account for the suppression of J/ip 
and T production in nuclear collisions in regions of phase space with high 
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particle densityEJ 

It is of particular interest to examine the fragmentation of the proton when 
the electron strikes a light quark and the interacting Fock component is the 
\uudcc) or \uudbb) state. These Fock components correspond to intrinsic charm 
or intrinsic bottom quarks in the proton wavcfunction. Since the heavy quarks 
in the proton bound state have roughly the same rapidity as the proton itself, 
the intrinsic heavy quarks will appear in the proton fragmentation region. 
One expects heavy quarkonium and also heavy hadrons to be formed from the 
coalescence of the heavy quark with the valence u and d quarks, since they 
have nearly the same rapidity. 



12 Intrinsic Charm and the J/ip — ► pit problem 

One of the most dramatic problems confronting the standard picture of quarko- 
nium decays is the J/iL^-> pn puzzleu This decay occurs with a branching 
ratio of (1.28 ± 0.10)%pl and it is the largest two-body hadronic branching 
ratio of the J/ip. The J/tp is assumed to be a cc bound state pair in the ^(lS 1 ) 
state. One then expects the tp' = ^(2S) to decay to pn with a comparable 
branching ratio, scaled by a factor ~ 0.15, due to the ratio of the ^(2S) to 
3L/1S 1 ) wavefunctions squared at the origin. In fact, B(ip' — > pn) < 3.6 x 10~ 5 
more than a factor of 50 below the expected rate. Most of the branching ra- 
tios for exclusive hadronic channels allowed in both J/ tp and ip' decays indeed 
scale with their lepton pair branching ratios, as would be expected from_decay 
amplitudes controlled by the quarkonium wavefunction near the originplLj 

B{%p' -> h) B(tb' ->■ e+e") , , 

; -0.147 ± 0.023 (66) 



B{J/ip^h) B(J/tp -> e+e~) 

where h denotes a given hadronic channel. The J/tp — > pn and J/tp — > KK* 
decays also conflict dramatically with perturbative QCD hadron helicity con- 
servation: all such pseudoscalar/vector two-body hadronic finaLstates are for- 
bidden at leading twist if helicity is conserved at each vertexEH 

Marek Kar liner and IE3 have recently shown that such anomalously large 
decay rates for the J/ip and their suppression for tp (2S) follow naturally from 
the existence of intrinsic charm | qq cc) Fock components of the light vector 
mesons. For example, consider the light-cone Fock representation of the p: 
p+ — \$ p _ I ud) + ^ p -_ I udcc) + ■ ■ ■ . The 4' p - wavefunction will be 

ud ' \ , ud cc [ ' udcc 

maximized at minimal invariant mass; i.e. at equal rapidity for the constituents 
and in the spin configuration where the ud are in a pseudoscalar state, thus 
minimizing the QCD spin-spin interaction. The cc in the | udcc) Fock state 
carries the spin projection of the p. We also expect the wavefunction of the cc 
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quarks to be in an S'-wave configuration with no nodes in its radial dependence, 
in order to minimize the kinetic energy of the charm quarks and thus also 
minimize the total invariant mass. 

The presence of the | udcc) Fock state in the p will allow the J/ip — » p-K 
decay to occur simply through rearrangement of the incoming and outgoing 
quark lines; in fact, the | udcc) Fock state wavefunction has a good overlap 
with the radial and spin | cc) and | ud) wavefunctions of the J/ \j) and pion. 
Moreover, there is no conflict with hadron helicity conservation, since the cc 
pair in the p is in the 1~ state. On the other hand, the overlap with the ip' will 
be suppressed, since the radial wavefunction of the n = 2 quarkonium state 
is orthogonal to the node-less cc in the | udcc) state of the p. This simple 
argument provides a compelling explanation of the absence of ip' — > pir and 
other vector pseudoscalar-scalar states]] 



13 Light-Cone Wavefunction Description of the Spin Anomaly in 
Deep Inelastic Polarized Structure Functions 

One of the most interesting distinguishing characteristics between extrinsic 
and intrinsic heavy quarks is their contributions to the Ellis-Jaffc sum rule 
/p 1 dxg\(x, Q) for polarized deep inelastic scattering cross sections. The ex- 
trinsic contributions to structure functions can be identified with photon-gluon 
fusion processes since they derive from QQ constituents of the gluon. However, 
one obtains zero contribution to the Ellis-Jaffe sum rule from j*g — * qq at tree 
level if the gluon is on-sheLUs 2 = 0. This follows from the DHG sum rule: the 
tree graph contribution tol£j 



/ — Aa(ab 



cd) (67) 



vanishes for any two-to-two polarized cross sections if a is an on-shell gauge 
particle. Thus the anomaly contribution AG to the Ellis-Jaffe sum rule 
arises from off-shell gluons with |fc 2 1 ^ 4mgL3 The final state which contributes 
physically to such configurations consists of Q and Q jets recoiling against the 
scattered lepton plus a third jet scattering at pr > rriQ, corresponding to a 
quark (or gluon) which emitted the off-shell gluon. The intrinsic contributions, 
on the other hand, consist of one high px heavy quark jet recoiling against the 



^ The possibility that the radial configurations of the initial and fmal states could be 
playing a role in the J/ip — > p-K puzzle was first suggested by S. PinskyO who however had 
in mind the radial wavefunctions of the light quarks in the p, rather than the wavefunction 
of the cc intrinsic charm components of the final state mesons. 
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lepton. Also, as noted above, the AQ(x) and AQ(x) in general are different 
hclicity distributionsta 

14 Direct Measurement of the Light-cone Valence Wavefunction. 

Diffractive multi-jet production in heavy nuclei provides a novel way to measure 
the shaoe-jOf the LC Fock state wavefunctions. For example, consider the 
reaction 

irA -> Jeti + Jet 2 + A' (68) 

at high energy where the nucleus A' is left intact in its ground state. The 
transverse momenta of the jets have to balance so that fcjj + k±2 = Q± < T^-a > 
and the light-cone longitudinal momentum fractions have to add to x\ + X2 ~ 1 
so that ApL < IZ^ 1 ■ The process can then occur coherently in the nucleus. 
Because of color transparency; i.e. the cancellation of color interactions in 
a small-size color-singlet hadron, the valence wavefunction of the pion with 
small impact separation,_anll penetrate the nucleus with minimal interactions, 
diffracting into jet pairsE3 The x\ = x, x^ = 1 — x dependence of the di-jet 
distributions will thus reflect the shape of the pion distribution amplitude; the 
fcj_l — k±2 relative transverse momenta of the jets also gives key information 
on the underlying shape of the valence pion wavefunction. The QCD analysis 
can be confirmed by the observation that the diffractive nuclear amplitude 
extrapolated to t — is linear in nuclear number A, as predicted by QCD color 
transparency. The integrated diffractive rate should scale as A 2 /1Z\ ~ A 4 / 3 . 
A diffractive experiment of this type is_now in progress at Fermilab using 500 
GeV incident pions on nuclear targets cB 

Data from CLEO for the 77* — > 7r° transition form factor favor a form for 
the pion distribution amplitude close to the asymptotic solution^ dtff vmpt (x) = 
\/3/xx(l — x) to the perturbative QCD evolution equationt3E2lO It will be 
interesting to see if the diffractive pion to di-jet experiment also favors the 
asymptotic form. 

It would also be interesting to study diffractive tri-jet production using pro- 
ton beams pA Jeti + Jet 2 + Jet3 + A' to determine the fundamental shape 
of the 3-quark structure of the valence light-cone wavefunction of the nucleon 
at small transverse separation. Conversely, one can use incident real and vir- 
tual photons: ^*A — > Jeti + Jet 2 + A' to confirm the shape of the calculable 
light-cone wavefunction for transversely-polarized and longitudinally-polarized 
virtual photons. Such experiments will open up a remarkable, direct window 
on the amplitude structure of hadrons at short distances. 
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15 Other Applications of Light-Cone Quantization to Hadron Phe- 
nomenology 

The light-cone formalism provides the theoretical framework which allows for 
a hadron to exist in various Fock configurations. For example, quarkonium 
states not only have valence QQ components but they also contain QQg and 
QQgg states in which the quark pair is in a color-octet configuration. Similarly, 
nuclear LC wave functions contain components in which the quarks are not in 
color-singlet nucleon sub-clusters. In some processes, such as large momentum 
transfer exclusive reactions, only the valence color-singlet Fock state of the 
scattering hadrons with small inter-quark impact separation b± = 0(l/Q) can 
couple to the hard scattering amplitude. In reactions in which large numbers of 
particles are produced, the higher Fock components of the LC wavefunction will 
be emphasized. The higher particle number Fock states of a hadron containing 
heavy quarks can be diffractively excited, leading to heavy hadron production 
in the high momentum fragmentation region of the projectile. In some cases the 
projectile's valence quarks can coalesce with quarks produced in the collision, 
producing unusual leading-particle correlations. Thus the multi-particle nature 
of the LC wavefunction can manifest itself in a number of novel ways. For 
example: 

Regge behavior. The light-cone wavefunctions ip n /H °f a hadron are not 
independent of each other, but rather are. coupled via the equations of mo- 
tion. Recently Antonuccio, Dalley and iEa have used the constraint of finite 
"mechanical" kinetic energy to derive "ladder relations" which interrelate the 
light-cone wavefunctions of states differing by 1 or 2 gluons. We then use these 
relations to derive the Regge behavior of both the polarized and unpolarized 
structure functions at x — > 0, extending Mueller's derivation of the BFKL hard 
QCD pomeron from the_properties of heavy quarkonium light-cone wavefunc- 
tions at large N c QCDEI 

Analysis of diffractive vector meson photoproduction. The light-cone Fock 
wavefunction representation of hadronic amplitudes allows a simple eikonal 
analysis of diffractive high energy processes, such as r j*(Q 2 )p — * pp, in terms 
of the virtual photon and the vector meson Fock state lightpcpne wavefunctions 
convoluted with the gp — ► gp near- forward matrix element O See Fig. lh. One 
can easily show that only small transverse size b± ~ 1 /Q of the vector meson 
wavefunction is involved. The hadronic interactions are minimal, and thus the 
7*(Q 2 )iV — > pN reaction can occur coherently throughout a nuclear target in 
reactions such as without absorption or shadowing. The 7* A — > 4>A process 
thus provides a natural framework for testing QCD color transparency^ 
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Structure functions at large Xbj ■ The behavior of structure functions where 
one quark has the entire momentum requires the knowledge of LC wavcfunc- 
tions with x — > 1 for the struck quark and x — > for the spectators. As 
mentioned in Section 2, this is a highly off-shell configuration, and thus one 
can rigorously derive quark-counting and helicity-retention rules for the power- 
law behavior of the polarized and unpolarized quark and gluon distributions 
in the x — > 1 endpoint domain. Evolution of structure functions is minimal 
in this domain because the struck quark is highly virtual as x — > 1; i.e. the 
starting point Qo f° r evolution caimpt be held fixed, but must be larger than 
a scale of order (to 2 + k^)/(l — x)E3 

Color Transparency QCD predicts that the Fock components of a hadron 
with a small /polor dipole moment can pass through nuclear matter without 
interactions I 66 !! 74 ! Thus in the case of large momentum transfer reactions, where 
only small-size valence Fock state configurations enter the hard scattering am- 
plitude, both the initial and final state interactions of the hadron states be- 
come negligible. Color Transparency can be measured though the nuclear 
dependence of totally diffractive vector meson production da/dt(j*A — > VA). 
For large photon virtualities (or for heavy vector quarkonium) , the small color 
dipole moment of the vector system implies minimal absorption. Thus, re- 
markably, QCD predicts that the forward amplitude j*A — > VA at t — > 
is nearly linear in A. One is also sensitive to corrections from the non- 
linear ^4-dependence of the nearly forward matrix element that couples two 
gluons to the nucleus, which is closely related to the nuclear dependence of 
the gluon structure function of the nucleuslI3 The integral of the diffractive 
cross section over the forward peak is thus predicted to scale approximately as 
A 2 /R\ ~ A 4 / 3 . Evidence for color transparency in quasi-elastic p leptoproduc- 
tion 7*4pj p°N(A— 1) has recently been reported by the E665 experiment at 
FermilablJj for both nuclear coherent and incoherent reactions. A test could 
also be carried out at very small t min at HERA, and would provide a striking- 
test of QCD in exclusive nuclear reactions. There is also evidence for QCD 
"color transparency" in quasi-elastic pp scattering in nucleillj In contrast to 
color transparency, Fock states with large-scale color, configurations interact 
strongly and with high particle number production. l3 

Hidden Color The deuteron form factor at high Q 2 is sensitive to wavefunc- 
tion configurations where all six quarks overlap within an impact separation 
b^i < 0(1/ Q); the leading power-law fall off predicted by QCD i^i^Q 2 ) = 
f(a s (Q 2 ))/(Q 2 ) 5 , where, asymptotically, f(a s (Q 2 )) cx a s (Q 2 ) 5+2 ^ElThe deriva- 
tion of the evolution equation for the deuteron distribution amplitude and its 
leading anomalous dimension 7 is given in Ref. E3 In general, the six-quark 
wavefunction of a deuteron is a mixture of five different color-singlet states. 
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The dominant color configuration at large distances corresponds to the usual 
proton-neutron bound state. However at small impact space separation, all 
five Fock color-singlet components eventually acquire equal weight, i.e., the 
deuteron wavefunction evolves to 80% "hidden color." The relatively large 
normalization of the deuteron form factor observed at large Q 2 points to siz- 
able hidden color contributions E2l 

Spin-Spin Correlations in Nucleon-Nucleon Scattering and the Charm 
Threshold One of the most striking anomalies in elastic proton^proton scat- 
tering is the large spin correlation A^^ observed at large anglesO At y/s ~ 5 
GeV, the rate for scattering with incident proton spins parallel and normal 
to the scattering plane is four times larger than that for scattering with anti- 
parallel polarization. This strong polarization correlation can be attributed 
to the onset of charm production in the intermediate state at this energycS 
The intermediate state \uuduudcc) has odd intrinsic parity and couples to the 
J = S = 1 initial state, thus strongly enhancing scattering when the incident 
projectile and target protons have their spins parallel and normal to the scat- 
tering plane. The charm threshold can also explain the anomalous change in 
color transparency observed at the same energy in quasi-elastic pp scattering. 
A crucial test is the observation of open charm production near threshold with 
a cross section of order of 1/xb. 

The QCD Van Der Waals Potential and Nuclear Bound Quarkonium The 
simplest manifestation of the nuclear force is the interaction between two heavy 
quarkonium states, such as the T(bb) and the J/ip(cc). Since there are no va- 
lence quarks in common, the dominant color-singlet interaction arises simply 
from the exchange of two or more gluons. In principle, one could measure 
the interactions of such systems by producing pairs of quarkonia in high en- 
ergy hadron collisions. The same fundamental QCD van der Waals potential 
also dominates the interactions of heavy quarkonia with ordinary hadrons and 
nuclei. The small size of the QQ bound state relative to the much larger 
hadron allows a systematic expansion of the gluonic potential using the op- 
erator product expansionEj The coupling of the scalar part of the interaction 
to large-size hadrons is rigorously normalized to the mass of the state via the 
trace anomaly. This scalar attractive potential dominates the interactions at 
low relative velocity. In this way one establishes that the nuclear force be- 
tween heavy quarkonia and ordinary nuclei is attractive and sufficiently strong 
to produce nuclear-bound quarkoniumEjo Recently, Miller and I have shown 
that the corrections to the gluon exchange potential from meson exchange con- 
tributions are relatively negligible, and we show_how deuteron targets can be 
used to measure the J/?/>-nucleon cross sectioned Navarra and I have shown 
that exclusive decays of B mesons at B factories such as the B — > J/ippA can 
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provide a sensitive search tool for finding possible J/?/;-baryon resonances 



16 Commensurate Scale Relations 

A critical problem in making reliable predictions in quantum chromodynamics 
is how to deal with the dependence of the truncated perturbative series on 
the choice of renormalization scale and scheme. For processes where only the 
leading and next-to-leading order predictions are known, the theoretical un- 
certainties from the choice of renormalization scale and scheme are often much 
larger than the experimental uncertainties. The uncertainties introduced by 
the conventions in the renormalization procedure are amplified in processes in- 
volving more than one physical scale such as jet observables and semi-inclusive 
reactions. In the case of jet production at e + e~ colliders, the jet fractions de- 
pend both on the total center of mass energy s and the jet resolution parameter 
y (which gives an upperbound ys to the invariant mass squared of each individ- 
ual jet), different scale-setting strategies can lead to very different behaviors 
for the renormalization scale in the small y region. In the case of QCD pre- 
dictions for exclusive processes such as the decay of heavy hadrons to specific 
channels and baryon form factors at large momentum transfer, the scale ambi- 
guities for the underlying quark-gluon subprocesses are even more acute since 
the coupling constant a s (/i) enters at a high power. Furthermore, since the 
external momenta entering an exclusive reaction are partitioned among the 
many propagators of the underlying hard-scattering amplitude, the physical 
scales that control these processes are inevitably much softer than the overall 
momentum transfer. 

The renormalization scale ambiguity problem can be resolved if one can 
optimize the choices ofjScale and scheme according to some sensible criteria. 
In the BLM procedure 0, the renormalization scales are chosen such that all 
vacuum polarization effects from the QCD (3 function are re-summed into the 
running couplings. The coefficients of the perturbative series are thus identi- 
cal to the perturbative coefficients of the corresponding conformally invariant 
theory with (3 = 0. The BLM method has the important advantage of "pre- 
summing" the large and strongly divergent terms in the PQCD series which 
grow as n\(a s f3o) n , i e tha. infrared renormalons associated with coupling con- 
stant renormalization. BeII Furthermore, the renormalization scales Q* in the 
BLM method are physicaLp. the sense that they reflect the mean virtuality of 
the gluon propagators. LJe!IE3'E3 In fact, in the av(Q) scheme, where the QCD 
coupling is defined from the heavy quark potential, the renormalization scale 
is by definition the momentum transfer caused by the gluon. 

A basic principle of renormalization theory is the requirement that re- 
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lations between physical observables must be independent of renormalization 
scale and scheme conventions to any fixed order of perturbation theoryo In this 
section, I shall discuss high precision perturbative predictions which have no 
scale or scheme ambiguities. These predictions, called "Commensurate Scale 
Relations," Ej are valid for any renormalizable quantum field theory and thus 
may provide a uniform perturbative analysis of the electroweak and strong 
sectors of the Standard Model. 

Commensurate scale relations relate observables to observables, and thus 
are independent of theoretical conventions such as choice of intermediate renor- 
malization scheme. The scales of the effective charges that appear in commen- 
surate scale relations are fixed by the requirement that the couplings sum all of 
the effects of the nonzero (3 function, as in the BLM methodd The coefficients 
in the perturbative expansions in the commensurate scale relations are thus 
identical to those of a corresponding conformally-invariant theory with /3 = 0. 

A helpful tool and notation for relating physical quantities is the effec- 
tive charge. Any perturbatpzabz calculable physical quantity can be used to 
define an effective charge lj'Ej'O by incorporating the entire radiative correc- 
tion into its definition. An important result is that all effective charges cka(Q) 
satisfy the Gell-Mann-Low renormalization group equation with the same (3q 
and (3\ ; different schemes or effective charges only differ through the third and 
higher coefficients of the (3 function. Thus, any effective charge can be used 
as a reference running coupling constant in QCD to define the renormalization 
procedure. More generally, each effective charge or renormalization scheme, 
including MS, is a special case of the universal coupling function a(Q, j3 n ). 
Peterman and Stiickelberg have shownE3 that all effective charges are related 
to each other through a set of evolution equations in the scheme parameters 

f3n. 

For example, consider the entire radiative corrections to the annihilation 
cross section expressed as the "effective charge" au(Q) where Q = s/s; 



R(Q) 



/ 



a R (Q) 



(69) 



Similarly, we can define the entire radiative correction to the Bjorken sum rule 
as the effective charge a gi (Q) where Q is the lepton momentum transfer: 



dxigr^Q^-gr^Q 2 )} 



2\1 — 



9A 




\ a gi (Q)' 


9v 




TT 



(70) 



The commensurate scale relations connecting the effective charges for ob- 
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servables A and B have the form 

cl a {Qa) = a B {Q B ) (l + r A/B — + •••) , (71) 

where the coefficient r^/B is independent of the number of flavors rip contribut- 
ing to coupling constant renormalization. We calculate the coefficients in the 
next section. The ratio of scales \a/b = Qa/Qb is unique at leading order and 
guarantees that the observables A and B pass through new quark thresholds 
at the same physical scale. One also can show that the commensurate scales 
satisfy the transitivity rule Xa/b — -^a/c-^c/Bj which is the renormalization 
group property which ensures that predictions in PQCD are independent of the 
choice of an intermediate renormalization scheme C. In particular, scale-fixed 
predictions can be made without reference to theoretically-constructed renor- 
malization schemes such as MSb. QCD can thus be tested in a new and precise 
way by checking that the observables track both in their relative normalization 
and in their commensurate scale dependence. 

A scale-fixed relation between any two physical observables A and B can 
be derived by applying BLM scale-fixing to their respective perturbative pre- 
dictions in, say, the MS scheme, and then algebraically eliminating ol-jj^. The 
choice of the BLM scale ensures that the resulting commensurate scale relation 
between A and, B is independent of the choice of the intermediate renormal- 
ization schemed Thus, using this formalism, one can relate any perturbatively 
calculable observable, such as the annihilation ratio R e + e - , the heavy quark 
potential, and the radiative corrections to structure function sup. rules to each 
other without any renormalization scale or scheme ambiguity^ Commensu- 
rate scale relations can also be applied in grand unified theories to make scale 
and scheme invariant predictions which relate physical observables in different 
sectors of the theory. 

Scales that appear in commensurate scale relations are physical since thpz. 
reflect the mean virtualityxif the gluons in the underlying hard subprocessQH 
As emphasized by Muellerlp commensurate scale relations isolate the effect of 
infrared renormalons associated with the nonzero /3 function. The usual facto- 
rial growth of the coefficients in perturbation theory due to quark and gluon 
vacuum polarization insertions is eliminated since such effects are resummed 
into the running couplings. The perturbative series is thus much more conver- 
gent. 

It is interesting to compare Pade resummation predictions for single-scale 
perturbative QCD series in which the initial renormalization scale choice is 
taken as tiue characteristic scale /i = Q as well as the BLM scale /i = Q*. 
One finds lj that the Pade predictions for the summed series are in each case 
independent of the initial scale choice, an indication that the Pade results are 
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thus characteristic of the actual QCD prediction. However, the BLM scale 
generally produces a faster convergence to the complete sum than the conven- 
tional scale choice. This can be understood by the fact that the BLM scale 
choice immediately sums into the coupling all repeated vacuum polarization 
insertions to all orders, thus eliminating the large (f3oa s ) n terms in the series 
as well, as the nl growth characteristic of the infrared renormalon structure of 

pqcdBb 



17 The Generalized Crewther Relation 

In 1972 Crewther derived a remarkable consequence of the operator product 
expansion for conformally-invariant gauge theory. Crewther's relation has the 
form 

35 = KB! (72) 

where S is the value of the anomaly controlling 7r° — > 77 decay, K is the value of 
the Bjorken sum rule in polarized deep inelastic scattering, and R' is the isovec- 
tor part of the annihilation cross section ratio er(e + e~ ^hadrons)/<7(e + e~ — » 
Since S is unaffected by QCD radiative correctionspl Crewther's rela- 
tion requires that the QCD radiative corrections to R e + e ~ exactly cancel the 
radiative corrections to the Bjorken sum rule order by order in perturbation 
theory. 

However, Crewther's relation is only valid in the case of conformally- 
invariant gauge theory, i.e. when the coupling a s is scale invariant. However, 
in reality the radiative corrections to the Bjorken sum rule and the annihila- 
tion ratio are in general functions of different physical scales. Thus Crewther's 
relation cannot be tested directly in QCD unless the effects of the nonzero (3 
function for the QCD running coupling are accounted for, and the energy scale 
y/s in the annihilation cross section is related to the momentum-transfer Q 
in the deep inelastic sum rules. Recently Broadhurst and KataevO have ex- 
plicitly calculated the radiative corrections to the Crewther relation and have 
demonstrated explicitly that the corrections are proportional to the QCD f3 
function. 



We can use the known expressions to three loops E£j£3E£j in MS scheme 
and choose the leading-order and next-to-leading scales Q* and Q** to re-sum 
all quark and gluon vacuum polarization corrections into the running couplings. 
The values of these scales are the physical values of the energies or momentum 
transfers which ensure that the radiative corrections to each observable passes 
through the heavy quark thresholds at their respective commensurate physi- 
cal scales. The final result connecting the effective charges (see Section 1) is 
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remarkably simple: 

a gi (Q) _ a R (Q*) fa a (Q**)Y , ( a R {Qf) 



I • (73) 

The coefficients in the series (aside for a factor of CV, which can be absorbed 
in the definition of a s ) are actually independent of color and are the same in 
Abelian, non Abelian, and conformal gauge theory. The non-Abelian structure 
of the theory is reflected in the scales Q* and Q** . Note that the a calculational 
device; it simply serves as an intermediary between observables and does not 
appear in the final relation ((73|). This is equivalent to the group property de- 
fined by Peterman and StiickelbergEHl which ensures that predictions in PQCD 
are independent of the choice of an intermediate renormalization sche me ( The 



renormalization group met hod was developed by Gell-Mann and LowL_3 and 



by Bogoliubov and ShirkovE_j) 



The connection between the effective charges of observables given by Eq. 
(f73|) is a prime example of a "commensurate scale relation" (CSR). A funda- 
mental test of QCD will be to verify empirically that the related observables 
track in both normalization and shape as given by the CSR. The commensu- 
rate scale relations thus provide fundamental tests of QCD which can be made 
increasingly precise and independent of the choice of renormalization scheme 
or other theoretical convention. More generally, the CSR between sets of phys- 
ica^ observables automatically satisfy the transitivity and symmetry properties 
of the scale transformations of the_renormalization "group" as originally 



10:, 



defined by Peterman and StiickelbergO The predicted relation between ob- 
servables must be independent of the order one makes substitutions; i.e. the 
algebraic path one takes to relate the observables. 

The relation-between scales in the CSR is consistent with the BLM scale- 
fixing procedure u in which the scale is chosen such that all terms arising from 
the QCD f3— function are resummed into the coupling. Note that this also im- 
plies that the coefficients in the perturbation CSR expansions are independent 
of the number of quark flavors / renormalizing the gluon propagators. This 
prescription ensures that, as in quantum electrodynamics, vacuum polarization 
contributions due to fermion pairs are all incorporated into the coupling a(/x) 
rather than the coefficients. The coefficients in the perturbative expansion 
using BLM scale-fixing are the same as those of the corresponding confor- 
mally invariant theory with /3 = 0. In practice, the conformal limit is defined 
by (3o,(3\ — > 0, and can be reached, for instance, by adding enough spin-half 
and scalar quarks as in N — 4 supersymmetric QCD. Since all the running 
coupling effects have been absorbed into the renormalization scales, the BLM 
scale-setting method correctly reproduces the perturbation theory coefficients 
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of the conformally invariant theory in the (3 — ► limit. 

The commensurate scale relation between a gi and o>r given by Eq. (73) 



implies that the radiative corrections to the annihilation cross section and the 
Bjorken (or Gross-Llewellyn Smith) sum rule cancel at their commensurate 
scales. The relations between the physical cross sections can be written in the 
forms: 



and 



3 E e 9 \9a/9v 

iW(s) J W 3 ^(x,Q 2 ) + jf (x,Q 2 ) 
3I> 2 6 



1 - Ap a 6 (74) 



= 1 - A/3 a 3 , (75) 



provided that the annihilation energy in R e + e - (s) and the momentum transfer 
Q appearing in the deep inelastic structure functions are commensurate at 
NLO: v/i = Q* = Qexp[| - 2C 3 + (& + fCa - Xl ~ £)A>a(Q)]. The light- 
by-light correction to the CSR for the Bjorken sum rule vanishes for three 
flavors. The term A/3 a 3 with A = in (Q** /Q*) is the third-order correction 
arising from the difference between Q** and Q*; in practice this correction is 
negligible: for a typical value a — Or(Q) /tt = 0.14, A/3 a 3 = 0.007. Thus at the 
magic energy ^/s = Q* , the radiative corrections to the Bjorken and GLLS sum 
rules almost precisely cancel the radiative corrections to the annihilation cross 
section. This allows a practical test and extension of the Crewther relation to 
nonconformal QCD. 

As an initial test of Eq. (|75|), we can compare the CCFR measurement 



of the Gross-Llewellyn Smith sum rule 1 — ap 3 — ^ f Q dx[F^ p (x,Q 2 
F^ p {x,Q 2 )] = |(^,± 0.13) at Q 2 = 3 GcV 2 and the parameterization of the 
annihilation datat23 1 + olr = R e + e - (s)/3 e q = 1-20. at the commensurate 
scale V* = Q* = 0.38 Q = 0.66 GeV. The product is (1 + a R )(l - S Fs ) = 
1.00 ±0.04, which is a highly nontdvial check of the theory at very low physical 



scales. More recently, the E143 £23 experiment at SLAC has reported a new 
value for the Bjorken sum rule at Q 2 = 3 GeV 2 : T%-F£ = 0.163±0.010(stat)± 
0.016(syst). The Crewther product in this case is also consistent with QCD: 
(l + a R )(l-a gi ) =0.93 ±0.11. 

In a paper with Gabadadze, Kataev and Lu □ we show that it is also 
possible and convenient to choose one unique mean scale Q in au(Q) so that 
the perturbative expansion will also reproduce the coefficients of the geometric 
progression. The possibility of using a single scale in the generalization of the 
BLM prescription beyond the next-to-leading order (NLO) was first considered 



by Grunberg and KataevtHa The new single-scale Crewther relation has the 
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form: 

a gi (Q) = a R (Q*) - a 2 R (Q*) + a 3 R (Q*) + ••-, (76) 

The generalized Crewther relation provides an important test of QCD. 
Since the Crewther formula written in the form of the CSR relates one ob- 
servable to another observable, the predictions are independent of theoretical 
conventions, such as the choice of renormalization scheme. It is clearly very 
interesting to test these fundamental self-consistency relations between the 
polarized Bjorken sum rule or the Gross-Llewellyn Smith sum rule and the 
e + e~-annihilation i?-ratio. Present data are consistent with the generalized 
Crewther relations, but measurements at higher precision and energies will be 
needed to decisively test these fundamental connections in QCD. 

It is worthwhile to point out that commensurate scale relations are derived 
within the framework of perturbation theory in leading twist and_dQ not involve 



the nonpcrturbative contributions to the Adler's function D(Q p3 and the R 



ratio, 
rules. 



112 



. 1 as to the polarized Bjorken and the Gross-Llewellyn Smith sum 
These nonperturbative contributions are expected to be significant 
at small energies and momentum transfer. In order to make these contributions 
comparatively negligible, one should choose relatively large values for s and 
Q 2 . In order to put the analysis of the experimental data for lower energies on 
more solid ground, it will be necessary to understand whether there exist any 
Crewther-type relations hetween nonperturbative order 0(l/Q 4 )-corrections 

In-order 0(1/ Q 2 ) higher twist contributions 



to the Adler's D-functionEL3 and 



to the deep-inelastic sum rules 

Commensurate scale relations such as the generalized Crewther relation 
discussed here open up additional possibilities for testing QCD. One can com- 
pare two observables by checking that their effective charges agree both in 
normalization and in their scale dependence. The ratio of leading-order com- 
mensurate scales Xa/b is fixed uniquely: it ensures that both observables A 
and B pass through heavy quark thresholds at precisely the same physical 
point. The same procedure can be applied to multi-scale problems; in gen- 
eral, the commensurate scales Q*, Q** , etc. will depend on all of the available 
scales. 

The coefficients in a CSR are identical to the coefficients in a conformal 
theory where explicit renormalon behavior does not appear. It is thus reason- 
able to expect that the series expansions appearing in the CSR are convergent 
when one relates finite observables to each other. Thus commensurate scale 
relations between observables allow tests of perturbative QCD with higher and 
higher precision as the perturbative expansion grows. 
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18 Renormalization Scale Fixing In Exclusive Processes 



As we have noted, perturbative QCD can be used to analyze a number of 
exclusive processes involving large momentum transfers, including the decay 
of heavy hadrons to specific channels such as B — > nir and T — > pp, baryon 
form factors at large t, andjised 9 c . m . hadronic scattering amplitudes such as 



7P — > Ti + n at high energieslLJj As in the-case of inclusive reactions, factoriza- 



Rff a 

tion theorems for exclusive processes 1131121 allow the analytic separation of the 
perturbatively-calculable short-distance contributions from the long-distance 
nonperturbative dynamics associated with hadronic binding. 

The scale ambiguities for the underlying quark-gluon subprocesses are par- 
ticularly acute in the case of QCD predictions for exclusive processes, since the 
running coupling a s enters at a high power. Furthermore, since each external 
momentum entering an exclusive reaction is partitioned among the many prop- 
agators of the underlying hard-scattering amplitude, the physical scales that 
control these processes are inevitably much softer than the overall momen- 
tum transfer. Exclusive process phenomenology is further complicated by the 
fact that the scales of the running couplings in the hard-scattering amplitude 
depend themselves on the shape of the hadronic wavefunctions. 

In this section we will discuss the application of the BLM method to fix 
the renormalization scale of the QCD coupling in exclusive hadronic ampli- 
tudes such as the pion form factor, the photon-to-pion transition form fac- 
tor and 77 — > 7r + 7r~ at large momentum transfer. Renormalization-schcme- 
independent commensurate scale relations will be established which connect 
the hard scattering subprocess amplitudes that control these exclusive pro- 
cesses to other QCD observables such as the heavy quark potential and the 
electron-positron annihilation cross section. Because the renormalization scale 
is small, we will argue that the effective coupling is nearly constanL-thus ac- 
counting for the nominal (dimensional counting) scaling behavior lii-3 of the 



data 



LI. 
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The heavy-quark potential V(Q 2 ) can be identified as the two-particle- 
irreducible scattering amplitude of test charges, i.e., the scattering of an in- 
finitely heavy quark and antiquark at momentum transfer t — —Q 2 . The rela- 
tion 

V(Q 2 ) = ~ 4nCF «: {Q2 \ (77) 

with Cf — {Nq — \)/2Nc — 4/3, then defines the effective charge ay{Q)- 
This coupling provides a physically-based alternative to the usual MS scheme. 
Recent lattice gauge calculations have provided strong constraints on the nor- 
malization and shape of ay(Q 2 ). 
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As in the corresponding case of Abelian QED, the scale Q of the coupling 
a v(Q) is identified with the exchanged momentum. All vacuum polarization 
corrections due to fermion pairs are incorporated in terms of the usual vacuum 
polarization kernels defined in terms of physical mass thresholds. The first 
two terms (3o = 11 — 2n//3 and (3\ = 102 — 38n//3 in the expansion of the f3 
function defined from the logarithmic derivative of ay(Q) are universal, i.e., 
identical for all effective charges at Q 2 ^> 4ra?. The coefficient /?2 for ay has 
recently been calculated in the MS schemetill 

The scale-fixed relation between ay and the conventional MS coupling is 

*ms(Q) = ay(e^Q) (l + ^4251 + • • •) , (78) 

above or below any quark mass threshold. The factor e 5 / 6 ~ 0.4346 is the ratio 
of commensurate scales between the two schemes to this order. It arises because 
of the convention used in defining the modified minimal subtraction scheme. 
The scale in the MS scheme is thus a factor ~ 0.4 smaller than the physical 
scale. The coefficient 2CU/3 in the NLO term is a feature of the non- Abelian 
couplings of QCD; the same coefficient occurs even if the theory had been 
conformally invariant with /3q = 0. The commensurate scale relation between 
ay, as defined from the heavy quark potential, and ajjg- provides an analytic 
extension of the MS scheme in which flavor thresholds are automatically, taken 



into account at their proper respective scalesliiS'Eil The coupling ay provides 
a natural scheme for computing exclusive amplitudes. Once we relate form 
factors to effective charges based on observables, there are no ambiguities due 
to scale or scheme conventions. 

The use of ay as the expansion parameter with BLM scale-fixing has 
also been found to be valuable in lattice gauge theory, greatly increasing the 
convergence of perturbative expansions relative to those using tha-bare lattice 



couplingEl In fact, new lattice calculations of the T spectruntLia have been 
used to determine the normalization of the static heavy quark potential and 
its effective charge: 

a^(8.2 GeV) = 0.196(3), (79) 
where the effective number of light flavors is rif — 3. The corresponding 



modified minimal subtraction coupling evolved to the Z mass using Eq. (78) 
is given by 

a Ws^ = °- 115 (2)- ( 8 °) 

This value is consistent with the world average of 0.117(5), but is significantly 
more precise. These results are valid up to NLO. 



48 



19 Hard Exclusive Two-Photon Reactions 



Exclusive two-photon processes such as 77 — > hadron pairs and the transi- 
tion form factor 7*7 — > neutral mesons play a unique role in testing quantum 
chromodynamics because of the simplicity of the initial stateE3 At large mo- 
mentum transfer the direct point-like coupling of the photon dominates at 
leading twist, leading to highly specific predictions which depend on the shape 
and normalization of the hadron distribution amplitudes 4>h{xit Q), the basic 
valence bound state wavefunctions. The most recent exclusive two-photon pro- 
cess data from CLEO provides stringent tests of these fundamental QCD 
predictions. 

Exclusive processes are particularly sensitive to the unknown nonperturba- 
tive bound state dynamics of the hadrons. However, in some important cases, 
the leading power-law behavior of an exclusive amplitude at large momentum 
transfer can be computed rigorously via a factorization theorem which sepa- 
rates the soft and hard dynamics. The key ingredient is the factorization of 
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aration of the perturbatively-calculable short-distance contributions from the 
long-distance nonperturbative dynamics associated with hadronic binding. For 
example, the amplitude 77 — * 7r + 7r~ factorizes in the form 

M ry ^ r +^- = I dx I dy (/>n(x,Q)T H (x,y,Q) (j)Tr(y,Q) (81) 
Jo Jo 

where ^(x, Q) is in the pion distribution amplitude and contains all of the soft, 
nonperturbative dynamics of the pion qq wavefunction integrated in relative 
transverse momentum up to the separation scale k\ < Q 2 , and Tjj is the 
quark/gluon hard scattering amplitude for 77 — > (qq)(qq) where the outgoing 
quarks are taken collinear with their respective pion parent. To lowest order 
in a s , the hard scattering amplitude is linear in a s . The most convenient 
definition of the coupling is the effective charge ay(Q 2 ), defined from the 
potential for the scattering of two infinitely heavy test charges, in analogy to 
the definition of the QED running coupling. Another possible choice is the 
effective charge aji(s), defined from the QCD correction to the annihilation 
cross section: R e + e - ^hadrons ( s ) = ^o(l + Q; i?( s )/ 7r )- One_can relate ay and 
an to cnjjg to NNLO using commensurate scale relationsEil 

The contributions from non-valence Fock states and the correction from 
neglecting the transverse momentum in the subprocess amplitude from the 
nonperturbative region are higher twist, i.e., power-law suppressed. The trans- 
verse momenta in the perturbative domain lead to the evolution of the distri- 
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bution amplitude and to next-to-leading-order (NLO) corrections in a s . The 
contribution from the endpoint regions of integration, x ~ 1 and y ~ 1, are 
power-law and SudakcpL-Suppressed and thus can only contribute corrections 
at higher order in 1/ QEJ 

The QCD coupling is typically evaluated at quite low scales in exclusive 
processes since the momentum transfers has to be divided among several con- 
stituents. In the BLM procedure, the scale of the coupling is evaluated by 
absorbing all vacuum polarization corrections with the scale of the coupling or 
by taking the experimental value integrating over the gluon virtuality. Thus, 
in the case of the (timelike) pion form factor the relevant scale is of order 
Q* 2 ~ e^ 3 A4 2 r7r ^ = t^j M 2 K+ ^- assuming the asymptotic form of the pion dis- 
tribution amplitude ^y m P t = \Z3f n x(l — x). At such low scales, it is likely 
that the coupling is frozen or relatively slow varying. 

In the BLM procedure, the renormalization scales are chosen such that all 
vacuum polarization effects from the QCD (3 function are re-summed into the 
running couplings. The coefficients of the perturbative series are thus identi- 
cal to the perturbative coefficients of the corresponding conformally invariant 
theory with (3 = 0. The BLM method has the important advantage of "pre- 
summing" the large and strongly divergent terms in the PQCD series which 
grow as n\(a s (3o) n , i.e_thc infrared renormalons associated with coupling con- 
stant renormalizationOO Furthermore, the renormalization scales Q* in the 
BLM method are physinaliii the sense that they reflect the mean virtuality of 
the gluon propagators EJ : I3'E30 In fact, in the av(Q) scheme, where the QCD 
coupling is defined from the heavy quark potential, the renormalization scale is 
by definition the momentum transfer caused by the gluon. Because the renor- 
malization scale is small in the exclusive 77 processes discussed here, we will 
argue that the effective coupling is ne arly constant, thus accounting for the 
nominal scaling behavior of the dataE£Ei3 

Ji, Pang, Robertson, and 10 have recently analyzed the pion transition 
form factor F 1 7 — ► ir° obtained from ej — > e'ir°, the timelike pion form 
obtained from e + e~ — > tt + tt, and the 77 — > it + tt~ processes, all at NLO 
in ay. The assumption of a nearly constant coupling in the hard scatter- 
ing amplitude at low scales provides an explanation for the phenomenolog- 
ical success of dimensional counting rules for exclusive processes; i.e., the 
power-law fall-off follows the nominal scaling of the hard scattering ampli- 
tude .Mhad ~ Th ~ [PpJn — 71 where n is in the total number of incident and 
final fields entering Th Eli 

The transition form factor has now been measured up to Q 2 < 8 GeV 2 
in the tagged two-photon collisions — > e'ir° by the CLEO and CELLO 
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collaborations. In this case the amplitude has the factorized form 



F lM (Q 2 ) = ^= J dxt M (x, Q 2 )T^ M (x, Q 2 ), (82) 
where the hard scattering amplitude for 77* — > qq is 

T 7 g M (x,Q 2 )= {1 _\ )Q2 (! + %))■ (83) 



The leading OCF) corrections have been computed by BraatenE 2 ^ and Dittes 
and Radyushkin ££3; however, the NLO corrections are necessary to fix the 
BLM scale at LO. Thus it is not yet possible to rigorously determine the BLM 
scale for this quantity. We shall here assume that this scale is the same as that 
occurring in the prediction for F^. For the asymptotic distribution amplitude 
we thus predict 

Q^(Q*) = 2/„(l-|^^). (84) 

As we shall see, given the phenomenological form of oty we employ (discussed 
below), this result is not terribly sensitive to the precise value of the scale. 

An important prediction resulting from the factorized form of these results 
is that the normalization of the ratio 



R *(Q 2 ) = JZFaTB (n 2\\2 ( 85 ) 



F V (Q 2 ) 
4wQ 2 \F W7 (Q 

= ajj^e-W-o) (1 _ 0.56^) (86) 

= a v {e- 3/2 Q) (l + 1-43^) (87) 

= aR (e 5 / 12 ^Q)(l-0M^) (88) 

is formally independent of the form of the pinn distribution amplitude. The 
a ~M~Q correction follows from combined referencesElillHj'H The next-to-leading 



correction given here assumes the asymptotic distribution amplitude. 

We emphasize that when we relate to ay we relate observable to ob- 
servable and thus there is no scheme ambiguity. Furthermore, effective charges 
such as oty are defined from physical observables and thus must be finite even 
at low momenta. A number of proposals have been suggested for the form of 
the QpQ coupling in the low-momentum regime. For example, Petronzio and 
Parisill2J have argued that the coupling must freeze at low momentum transfer 
in order that perturbative QCD loop integrations be well defined. Mattingly 
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and StevcnsonE£ZI have incorporated such behavior into their parameterizations 
of an at low scales. GribovEHj has presented novel dynamical arguments ^& 



lated to the nature of confinement for a fixed coupling at low scales. Zerwafii^ 



has noted the heavy quark potential must saturate to a Yukawa form since 
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at large distances. Cornwall E£D and others have argued that the gluon 
propagator will acquire an effective gluon mass m g from nonperturbative dy- 
namics, which again will regulate the form of the effective couplings at low 
momentum. We shall adopt the simple parameterization 



a v (Q) 



in 



(89) 



which effectively freezes the ay effective charge to a finite value for O 2 < 4m 2 ,. 

We can use the nonrelativistic heavy quark lattice results Eijt£3 to fix 
the parameters. A fit to the lattice data of the above parameterization gives 
Ay = 0.16 GeV if we use the well-known momentum-dependent n/E£j Further- 



more, the value m 2 = 0.19 GeV 2 gives consistency with the frozen value of an 
advocated by Mattingly and Stevenson. 1123 Their parameterization implies the 
approximate constraint an{Q)/iT ~ 0.27 for Q = y/s < 0.3 GeV, which leads 
to ay (0.5 GeV) ~ 0.37 using the NLO commensurate scale relation between 
ay and aR. The resulting form for ay is shown in Fig. The corresponding 
predictions for an and a-jjg using the CSRs at NLO are also shown. Note 
that for low Q 2 the couplings, although frozen, are large. Thus the NLO and 
higher-order terms in the CSRs are large, and inverting them perturbatively 
to NLO does not give accurate results at low scales. In addition, higher-twist 
contributions to ay and a^, which are not reflected in the CSR relating them, 
may be expected to be important for low Q 2 



J.:; 2 



It is clear that exclusive processes such as the photon-to-pion transition 
form factors can provide a valuable window for determining the magnitude 
and the shape of the effective charges at quite low momentum transfers. In 
particular, we can check consistency with the ay prediction from lattice gauge 
theory. A complimentary method for determining ay at low momentum isj-tq 
use the angular anisotropy of e + e~ — > QQ at the heavy quark thresholds 
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It should be emphasized that the parameterization (|8^) is just an approximate 
form. The actual behavior of ay(Q 2 ) at low Q 2 is one of the key uncertainties 
in QCD phenomenology 

As we have emphasized, exclusive processes are sensitive to the magnitude 
and shape of the QCD couplings at quite low momentum transfer: Qy 2 ~ 
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e- 3 Q 2 ~ Q 2 /20 and Q^ 2 ~ Q 2 /5oEii The fact that the data for exclusive 
processes such as form factors, two photon processes such as 77 — * tt + tt~ , and 
photoproduction at fixed 9 c .m. are consistent with the nominal scaling of the 
leading-twist QCD predictions (dimensional counting) at momentum transfers 
Q up to the order of a few GeV can be immediately understood if the effective 
charges ay and an are slowly varying at low momentum. The scaling of 
the exclusive amplitude then follows that of the subpiocess amplitude Th with 
effectively fixed coupling. Donnachie and LandshoffE2j have also argued that a 
frozen coupling is needed to explain the observed t~ s scaling of da/dt {pp —> pp) 
at large s 3> — t. Note also that the Sudakov effect of the end point region 
is the exponential of a double log series if the coupling is frozen, and thus is 
strong. 
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Figure 6: The 7 — > n° trans itio n form factor. The solid line is the full prediction including 
the QCD correction [Eq. (BQ)]; the dotted line is the LO prediction Q 2 F 1 - K (Q 2 ) = 2/ w . 



In Fig. ||, we compare the recent CLEO datatH2l for the photon to pion 
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transition form factor with the prediction 



J F 77r (Q 2 ) =2/^1 



5a v (e-^ 2 Q) 

3 TT 



(90) 



The flat scaling of the Q 2 F y7I (Q 2 ) data from Q 2 = 2 to Q 2 = 8 GeV 2 pro- 
vides an important confirmation of the applicability of leading twist QCD to 
this process. The magnitude of Q 2 F 1 t T (Q 2 ) is remarkably consistent with the 
predicted form, assuming the asymptotic distribution amplitude and including, 
the III QCD radiative correction with a v (e~ 3/2 Q)/7r ~ 0.12. Radyushkin^ 
OngEHl and KrollEJ have also noted that the scaling and normalization of the 



photon-to-pion transition form factor tends to favor the asymptotic form for 
the pion distribution amplitude and rules out broader, distributions such as 
the two-humped form suggested by QCD sum rules. One cannot obtain a 
unique solution for the nonperturbative wavefunction from the F^ 1 data alone. 
However, we have the constraint that 



1 



1 5 a v (Q*) 

3 7T 



0.8 



(91) 



(assuming the renormalization scale we have chosen in Eq. (84) is approxi- 
mately correct). Thus one could allow for some broadening of the distribution 
amplitude with a corresponding increase in the value of ay at low scales. 

We have also analyzed the 77 — > ir + ir~ , K + K~ data. These data exhibit 
true leading-twist scaling (Fig. fj]) , so that one would expect this process to be 
a good test of theory. One can show that to LO 



4|fW(«) 



f (77- 



H+fi ) 



(92) 



in the CMS, where dt = (s/2)d(cos9 c , m .) and here F % (s) is the time-like pion 
form factor. The ratio of the time- like to space- like pion form factor for the 
asymptotic distribution amplitude is given by 



|F ftimclikc) ( _ Q2)| 
^(spacelike)/Q 2 A 



K(-Q* 2 )| 
av(Q* 2 ) 



(93) 



If we simply continue Eq. (|89| ) to negative values of Q 2 then for 1 < Q 2 < 10 
GeV 2 , and hence 0.05 < Q* 2 < 0.5 GeV 2 , the ratio of couplings in Eq. S 
is of order 1.5. Of course this assumes the analytic application of Eq. (^9|). 
Thus if we assume the asymptotic form for the distribution amplitude, then 
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wc predict 



(timclikc) 



(-Q 2 ) ~ (0.3 GeV 2 )/Q 2 and hence 



1(77 



.36 



1 



f (77^M+M") 



1 - cos 4 B r 



The resulting prediction for the combined cross section (7(77 



(94) 



K+K-% 



is shown in Fig. |?j, along with CLEO datatH3 Considering the possible contri 
bution of the resonance /2(1270), the agreement is reasonable. 




0.01 1 1 1 1 l 

12 3 4 

M (GeV/c 2 ) 



Figure 7: Two-photon annihilation cross section (7(77 — > 7r+7r , K + K ) as a function of 

CMS energy, for | cos6»*| < 0.6. 

We also note that the normalization of ay could be larger at low mo- 
mentum than our estimate. This would also imply a broadening of the pion 
distribution amplitude compared to its asymptotic form since one needs to 
raise the expectation value of 1/(1 — x) in order to maintain consistency with 
the magnitude of the F 1 - K (Q 2 ) data. A full analysis will then also require 

II The contribution from kaons is obtained at this order simply by rescaling the prediction 
for pions by a factor (/k//tt) 4 — 2.2. 
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consideration of the breaking of scaling from the evolution of the distribution 
amplitude. In any case, we find no compelling argument for significant higher- 
twist contributions in the few GeV regime from the hard scattering amplitude 
or the endpoint regions, since such corrections violate the observed scaling 
behavior of the data. 

The analysis we have presented here suggests a systematic program for 
estimating exclusive amplitudes in QCD (including exclusive B-decays) which 
involve hard scattering. The central input is ay(0), or 

Q 2 

aV=-L I ° dQ /2 a v {Q' 2 ), Q% < 1 GeV 2 , (95) 
Qo Jo 

which largely controls the magnitude of the underlying quark-gluon subpro- 
cesses for hard processes in the few-GeV region. In this work, the mean cou- 
pling value for Q\ ~ 0.5 GeV 2 is ay ~ 0.38. The main focus will then be to 
determine the shapes and normalization of the process-independent meson and 
baryon distribution amplitudes. 

The leading-twist scaling of the observed cross sections for exclusive two- 
photon processes and other fixed 9 cm reactions can be understood if the effec- 
tive coupling ay(Q*) is approximately constant in the domain of Q* relevant 
to the underlying hard scattering amplitudes. In addition, the Sudakov sup- 
pression of the long-distance contributions is strengthened if the coupling is 
frozen because of the exponentiation of a double log series. We have also found 
that the commensurate scale relation connecting the heavy quark potential, as 
determined from lattice gauge theory, to the photon-to-pion transition form 
factor is in excellent agreement with — ► ir°e data assuming that the pion 
distribution amplitude is close to its asymptotic form \Z3fnx(l— x). We also re- 
produce the scaling and approximate normalization of the 77 — ► 7r + 7r~, K + K~ 
data at large momentum transfer. However, the normalization of the space- 
like pion form factor F 7T (Q 2 ) obtained from electroproduction experiments is 
somewhat higher than that predicted by the corresponding commensurate scale 
relation. This discrepancy may be due to systematic errors introduced by the 
extrapolation of the j*p — > 7r + n electroproduction data to the pion pole. 
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